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ABSTRACT 

This  work  reviews  the  existing  literature  on  slam  force 
prediction  methods  and  discusses  some  of  the  theory  of  ship 
bottom  slamming.   One  of  these  prediction  methods  is  chosen 
as  the  best  and  further  improvements  are  made  on  it  by  using 
a  new  model  of  slamming. 

The  work  is  designed  to  determine  a  slam  forcing  function 
suitable  for  studying  ship  hull  girder  dynamic  response. 
Slam  forces  are  determined  by  integrating  the  pressure  time 
history  over  the  portion  of  the  hull  subject  to  slamming. 
A  computer  program  is  included  which  determines  the  slam 
forces  and  impulses  from  the  ship's  lines,  principal  dimen- 
sions, and  seakeeping  data.   The  method  is  applied  to  two 
ships,  the  MARINER  and  the  FOTINI-L,  and  the  results 
presented. 
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I.   INTRODUCTION 

The  primary  objective  of  this  thesis  is  to  review  the 
existing  literature  on  slamming  and  to  develop  a  method  of 
predicting  slam  forces  to  determine  hull  girder  response. 
This  project  is  a  portion  of  an  overall  study  presently  under- 
way at  M.I.T.  sponsored  by  the  American  Bureau  of  Shipping 
(A.B.S.)  on  hull  stiffness  criteria. 

Presently,  explicit  criteria  for  ship  stiffness  do  not 
exist.   The  ship  classification  societies  have  controlled  stiff- 
ness or  deflections  by  limiting  allowable  Length  to  Depth  (L/D) 
ratios  [1]*.   With  the  introduction  of  new  materials  such  as 
high  strength  steels  and  aluminum  in  ship  construction,  the 
question  of  stiffness  and  what  should  be  allowed  becomes  more 
important. 

Considering  the  present  L/D  ratios  as  hull  girder  stiff- 
ness criteria,  what  other  factors  that  affect  stiffness  could 
be  used  in  design?  Bending  stiffness  of  a  ship's  hull  can  be 
varied  within  reasonable  limits  by  the  following  means: 

1.  By  changing  the  basic  hull  depth, 

2.  By  designing  to  higher  stress  levels  with  higher 
strength  steels,  and 

3.  By  using  materials  with  different  elastic  moduli. 
Since  other  factors  can  more  properly  be  used  to  control  stiff- 
ness other  than  L/D,  it  is  not  satisfactory  to  use  this  as  a 


♦Numbers  in  square  brackets  designate  References  at  the  end 
of  the  paper. 
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criterion  for  design.   Probably  a  new  criterion  based 
specifically  on  a  measure  of  stiffness  should  be  used. 

Before  any  quantitative  information  could  be  developed 
for  limiting  hull  stiffness,  factors  which  would  be  changed 
by  stiffness  variations  were  assessed  [2] .   Assuming  strength 
criteria  are  maintained,  stiffness  variations  will  result  in 
changing  factors  which  determine  hull  girder  dynamic  response. 
In  the  Reference  [2]  report  it  was  concluded  that  the  factors 
of  ship  stiffness  which  should  be  of  concern  were: 

1.  Slamming  response  for  its  contribution  of  stress 
components,  especially  amidships, 

2.  springing  for  its  cumulative  stress  effects,  and 

3.  steady-state  propeller-excited  vibratory  motions 
and  their  deleterious  effects  upon  personnel  and 
main  machinery  components. 

Although  propeller-excited  vibratory  motions  will  change  with 
hull  stiffness  variations,  it  is  felt  that  these  problems  can 
be  dealt  with  by  more  efficient  and  less  costly  methods  than 
changing  the  hull  girder  stiffness.   Therefore  the  primary 
factors  of  concern  are  slamming  and  springing  response  with 
variations  in  hull  stiffness. 

In  the  overall  study  the  Kline-Clough  lumped-mass  vibra- 
tion program,  S.H.V.R.S.,  has  been  used  to  analyze  the  ship 
response  to  slamming  and  springing  [3] .   Several  ships  have 
been  looked  at  in  great  detail.   Some  of  the  results  have  been 
reported  in  references  [3]  and  [4].   To  date,  in  this  work  on 
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slamming,  a  constant  slam  impulse  of  100  Ton-seconds  has  been 
used  in  the  response  analysis  modeled  as  a  half  sine  pulse 
applied  at  a  single  location.   Results,  such  as  stress  levels, 
other  than  trends  in  response  due  to  changes  in  stiffness, 
could  not  be  determined  until  a  better  model  of  the  slam  force 
and  a  satisfactory  means  of  evaluating  force  magnitudes  from 
ship  to  ship  could  be  found. 

The  purpose  of  this  work  is  to  develop  a  model  of  the 
bottom  impact  slam  force  compatible  with  the  S.H.V.R.S.  program 
which  will  require  a  minimum  of  manhours  and  computer  time  yet 
give  adequate  results.   The  bottom  impact  slam  force  developed 
here  is  a  function  of  hull  form,  ship  speed,  and  sea  state. 
A  computer  program  is  developed  which  predicts  the  most 
probable  slam  force  as  well  as  an  extreme  force  to  be  used  for 
design.   Other  slam  information  is  also  available  such  as 
frequency  of  slamming  and  slam  pressures  useful  in  designing 
hull  plating. 
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II.   BOTTOM  IMPACT  SLAMMING 

For  this  study  the  bottom  impact  slam  is  defined  as  the 
"heavy  blow"  resulting  from  the  forebody  hitting  the  water 
surface  at  a  high  relative  velocity.   This  is  in  contrast  to 
the  phenomenon  of  bow  flare  slamming  which  is  a  consequence 
of  emersion  of  the  bow  flare  and  the  resultant  hydrodynamic 
and  hydrostatic  forces.   Much  of  the  work  to  date  has  been  in 
the  form  of  predicting  the  local  pressures  developed  to  study 
plating  damage.   That  can  be  considered  the  micro-problem.   In 
this  work  we  are  concerned  with  the  macro-problem  which  in- 
volves hull  vibrations  and  possibly  large  midship  stresses. 

Oakley  [5]  has  analyzed  the  micro-problem  in  great  detail. 
There  are  many  factors  that  have  been  observed  to  measurably 
affect  the  pressure  history.   These  factors  include  the 
following:   shape,  velocity,  and  phasing  of  the  surface  waves; 
geometry  of  the  hull  and  hull  reaction  characteristics;  air 
and  water  compressibility;  two  and  three-dimensional  air  and 
water  flow;  bubbles,  spray,  and  even  such  things  as  marine 
growth.   Due  to  the  complexity  of  the  problem  and  the  random- 
ness of  the  sea,  model  testing  and  statistics  play  an  important 
role  in  the  solution  to  the  slamming  problem. 

Tick  [6]  developed  the  equation  for  the  frequency  of 
slamming  which  is  dependent  on  the  following: 

1.  forefoot  emergence  statistics, 

2.  relative  velocity  statistics,  and 

3.  the  angle  between  the  keel  line  and  the  wave 
surface  at  the  instant  of  impact. 
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In  most  work  to  date  the  equation  is  simplified  to  remove  the 
angle  dependence  and  only  the  requirements  of  (1)  forefoot 
emergence  and  (2)  a  relative  velocity  greater  than  a  threshold 
velocity  are  considered  necessary  for  slamming.   Neglecting 
this  angular  dependence   esults  in  overpredicting  the  frequency 
of  slam  impact  as  well  a^  the  slam  force.   A  method  of  dealing 
with  the  angular  dependence  is  presented  in  Section  IV  of  this 
work.   The  two  requirements  for  slamming  are  met  by  the  ship 
in  Figure  1.   Using  Tick's  equation  without  angular  dependence 
allows  prediction  of  the  frequency  of  slamming  and  is  the  basis 
of  slam  calculations  in  seakeeping  programs  [7] .   The  one 
necessary  input  to  seakeeping  programs  is  the  threshold  velo- 
city.  How  this  is  detei mined  and  the  value  to  be  used  is 
discussed  in  Section  IV. 

From  Tick's  work  it  is  possible  to  conclude  whether  or 
not  a  ship  is  subject  to  slamming  under  a  particular  set  of 
operating  conditions,  in  particular  the  ship  speed,  sea  state 
and  draft  at  the  station  under  consideration.   Given  that 
slamming  does  occur,  the  next  step  is  to  determine  the  slam 
force.   Two  basic  methods  have  been  used  in  the  past:   (1) 
integration  of  the  slam  pressure  time  histories  over  the  slam 
area;  and  (2)  determine  an  analytical  solution  based  on  the 
rate  of  change  of  momentum.   In  Section  III  both  methods  are 
discussed  and  one  selected  for  this  study. 


15 


PHASE    (I) 


»HASE    (ll) 


PHASE     (III) 


FIGURE     2    -    PHASES     OF     SLAMMING 
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Before  comparing  the  two  methods,  some  discussion  is  in 
order  with  respect  to  the  mechanism  of  slamming.   Oakley  [5] 
considered  the  time  history  of  slamming  to  be  made  up  of 
three  phases: 

1.  Phase  (I)  is  the  time  in  which  the  body  is 
approaching  the  free  surface  up  until  it  makes 
contact.   In  this  phase  the  air  flow  and  the 
surface  wave  deflection  are  of  predominant 
importance . 

2.  Phase  (II)  is  the  period  when  the  body  goes 
through  the  free  surface  and  involves  the  problem 
of  "contact"  and  more  or  less  complete  wetting  of 
the  hull. 

3.  Phase  (III)  is  the  fully  wetted  problem  in  which 
the  rate  of  change  of  added  mass  and  the  local 
surface  deflection  are  dominant. 

A  very  flat  body  or  "U"  form  hull  may  accelerate  the  free 
surface  significantly  during  phase  (I)  while  a  more  "V" 
shaped  section  may  develop  greatest  forces  only  during  phase 
(III).   Figure  2  shows  the  relationship  between  the  three 
phases.   Much  of  the  work  to  date  has  dealt  with  only  one  or 
two  of  these  phases  at  a  time.   This  is  especially  true  with 
the  many  two  and  three-dimensional  drop  tests  performed. 
Chuang  [8]  and  others  have  demonstrated  the  important  role 
of  the  air  in  accelerating  the  free  surface  in  phase  (I)  and 
as  a  cushioning  mechanism  in  phase  (II) .   The  wedge-entry 
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problem  has  been  studied  by  many  authors,  most  recently  by 
Chuang  [9] ,  and  agreement  for  large  deadrise  angles  (phase 
(III)  only)  is  good.   The  theory  breaks  down  for  low  deadrise 
angles  where  phase  (I)  and  phase  (II)  become  important. 

Therefore  any  solution  to  the  slam  forcing  function 
should  take  into  account  the  forces  developed  during  all 
three  of  the  phases  unless  one  phase  can  be  shown  to  be  much 
greater  than  the  others.   At  this  time  it  appears  that  neglect- 
ing phase  (I)  and  phase   II)  is  only  valid  for  "V"  form  sections 
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III.   SLAM  FORCE  PREDICTION  METHODS 

During  the  last  twenty  years  a  great  deal  has  been  written 
on  slamming,  attesting  to  the  difficulty  of  a  satisfactory 
solution  to  the  problem.   Much  of  the  work  has  been  done  along 
with  ship  vibration  anal'sis  techniques.   This  study  does  not 
go  into  the  question  of  structural  response.   It  is  felt  that 
adequate  methods  exist  for  analyzing  the  structural  response 
of  a  ship  if  the  various  forcing  functions  can  be  adequately 
represented.   The  S.H.V.R.S.  program  described  in  Reference  [3] 
which  has  been  used  in  this  overall  study  has  been  shown  to 
give  satisfactory  results.   Therefore  the  model  of  slam  force 
was  chosen  with  the  S.H.V.R.S.  program  in  mind. 

The  slam  force  prediction  methods  are  of  two  types: 

1.  The  slam  force  is  considered  as  resulting  from 
the  rate  of  change  of  momentum  of  the  added 
mass  of  the  ship  with  a  buoyancy  correction 
term.   This  method  will  be  called  the  momentum 
method. 

2.  The  second  method  consists  of  determining  the  slam 
pressure  and  integrating  spatially  as  well  as 
temporally  over  the  slam  area.   This  method  will 
be  designated  the  slam  pressure  method. 

Both  of  these  methods  were  considered  and  the  second  was 
considered  more  valid  and  more  easily  applicable  to  the 
overall  stiffness  study. 
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RATE  OF  CHANGE  OF  MOMENTL  1  METHOD 

As  early  as  1956  Szebehely  [10]  proposed  that  the  slam 
force  was  a  result  of  the  rate  of  change  of  momentum  of  the 
added  mass  in  the  region  of  the  slam.   Further  studies  by 
Leibowitz  [11]  used  this  model  of  slamming  to  analyze  the 
slam  forces  on  a  Dutch  destroyer.   Good  results  were  obtained 
between  the  theory  and  full  scale  results  for  bow  flare  slams. 
More  recently  Mansour  anc  d'Oliviera  [12],  [13]  have  utilized 
the  same  slam  force  function  to  predict  slam  bending  moments 
in  regular  waves. 

The  slam  force  (P)  is  considered  as  the  sum  of  an  inertial 
term  (P.)  and  a  buoyancy  term  (P_)  and  is  given  by 


P(x,t)  =  P1(x,t)  +  P2(x,t) 


§¥(mwr)  +  PgA  (1) 


Tt(mWr}  +  lx(m"r)  aT  +  PgA 


where  m  =  m(x,t)  is  the  added  mass  per  unit  length,  w   is  the 
relative  vertical  velocity  between  the  wave  and  the  ship,  and 
A  is  the  sectional  area.   The  total  derivative  is  used  (— ) 
since  the  rate  of  change  of  momentum  occurs  not  only  temporally 
but  also  spatially. 

It  should  be  noted  that  this  is  the  phase  (III)  problem 
described  by  Oakley.  No  account  is  made  for  the  force  which 
comes  during  phase  (I)  and  phase  (II)  of  the  slam  process.   It 
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is  expected  that  this  method  should  show  good  results  for  "V" 
form  hulls  and  bow  flare  slamming  only.   The  major  drawback 
to  this  approach,  even  if  it  were  valid  for  the  bottom  impact 
slam  case,  is  it  requires  a  large  amount  of  seakeeping  infor- 
mation to  determine  m(x,t).   In  reality  the  added  mass  must 
be  determined  <is  a  function  of  draft  at  each  station  and 
frequency  of  encounter  f c r  each  wave  considered.   To  obtain 
this  information  requires  a  large  expenditure  in  manhours  and 
computer  time.   The  end  lesult  is  also  valid  only  for  regular 
waves . 

SLAM  PRESSURE  PREDICTION  METHOD 

A  great  deal  of  the  initial  work  in  this  area  has  been 
done  to  solve  the  micro-problem — that  is,  to  determine  the 
slam  pressure  locally  foi  plate  design  information.   Much  of 
the  work  in  this  area  ha:  been  done  by  M.D.  Ochi,  K.M.  Ochi, 
and  S.L.  Chuang  at  the  NU  val  Ship  Research  and  Development 
Center.   A  summary  of  the  slamming  studies  to  date  appears  in 
Reference  [14].   In  this  work  by  M.K.  Ochi  and  Motter  they 
develop  the  basic  theory  behind  the  approach  to  slamming  pre- 
diction in  this  method. 

In  this  approach  the  pressure  is  considered  as  a  function 
of  velocity  of  impact  and  section  shape  only.   The  slam 
pressure  is  given  by 


p  =  kf2  (2) 
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where  k  is  a  parameter  to  be  determined  from  section  form 
and  r  is  the  velocity  of  Impact. 

From  the  many  two  and  three-dimensional  drop  tests  as 
well  as  full  scale  measur sments  and  model  tests  it  has  been 
found  that  equation  (2)  h  )lds  quite  well  for  slam  impact  pro- 
blems.  Therefore  once  th 2  form  coefficient  (k)  can  be  deter- 
mined for  each  section,  tie  slam  impact  pressure  can  be 
determined  if  the  velocity  of  impact  is  known.   Since  relative 
velocities  are  easily  determined  from  seakeeping  programs, 
knowing  the  value  of  k  for  each  section  allows  for  the  deter- 
mination of  the  slam  pressure  which  can  be  integrated  over  the 
hull  form  to  determine  the  slam  force. 

This  second  approach  has  the  advantage  that  the  form 
coefficient  k  is  obtained  from  model  data  which  includes  the 
contribution  from  phase  (I),  (II),  and  (III)  forces.   In 
addition,  recent  work  by  Loukakis  and  Chryssostomidis  [15] 
make  available  the  seakeeping  information  necessary  for  slam 
studies  for  ships  that  can  be  approximated  by  Series  60  hull 
forms . 

The  slam  pressure  approach  was  chosen  since  it  requires 
a  minimum  of  seakeeping  information,  is  applicable  to  all 
hull  forms,  and  the  results  are  applicable  in  irregular  seas. 
Unlike  most  engineering  problems  it  also  turns  out  that  this 
approach  is  easiest  to  apply. 
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IV.   PRINCIPLES  AND  PROCEDURE  OF  THE  SLAM  PRESSURE  METHOD 

Ochi  and  Motter ' s  work  described  in  Reference  [14]  pre- 
sents the  theory  and  assumptions  of  the  slam  pressure  method. 
In  trying  to  apply  the  method  it  appeared  that  some  of  their 
assumptions  resulted  in  overpredicting  the  slam  force  and  im- 
pulse.  In  this  section  the  theory  and  assumptions  made  by 
Ochi  and  Motter  are  discussed  and  modified  where  the  assump- 
tions were  not  felt  full^  justified.   The  following  topics 
will  be  discussed  in  developing  the  method: 

1.  Frequency  of  slci  impact, 

2.  Slam  pressure  ar d  hull  form  coefficient, 

3.  Most  probable  extreme  pressure, 

4.  Extreme  pressure  for  design, 

5.  Sea  state  data  i  nd  approximation  for  storm  duration, 

6.  Slam  impact  are£  of  the  hull,  and 

7.  Slam  impact  force  and  impulse. 

FREQUENCY  OF  SLAM  IMPACT 

Tick  [6]  derived  the  equation  for  the  frequency  of  slam 
impact  at  the  bow  which  is  dependent  on  the  relative  velocity, 
relative  motion,  and  angle  between  the  keel  and  the  wave  sur- 
face.  By  integrating  over  all  angles  to  neglect  the  angular 
dependence, 


/r: 

N      =   tt-       f  ttt  ?r      {slam   impact} 

s        2ir      /     w  r 


„,2  .2 

(—  +  — )                                                  (3) 

(R'  R'.}                                                      K     ' 

r  r 
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where    N  =  Number  of  slams  per  second 
s 

T  =  Ship  draft  at  station  under  consideration 
r^  =  Threshold  velocity 
R'  =  Twice  variance  of  relative  motion 
R!  =  Twice  variance  of  relative  velocity 
The  number  of  impacts  in  a  storm  of  duration  t  hours,  N(t), 
is  then  given  by 


3600  t    /  "f         rT2    rj*,  .„, 

N(t)  =  — 5^ —   /  rT   exp  -[^-r  +  -f]  (4) 

r    r 


Different  values  of  the  threshold  velocity  below  which 
slamming  does  not  occur  have  been  used,  but  following  Ochi * s 
work  [14], 


**  =  12/5l0   =  °*53   ^  (5) 


will  be  used  where  r^  =  12  feet/sec  for  a  520  FT  ship  and  the 
Froude  scaling  law  is  applied.   Because  some  of  the  slams 
predicted  by  equation  (4)  will  occur  above  the  angle  which  will 
cause  slamming,  equation  (4)  will  overpredict  the  number  of 
slams.   This  number  is  used  to  determine  the  extent  of  slam 
area  and  slam  pressure.   Overestimating  N(t)  will  yield  a 
larger  slam  area  as  well  as  a  slightly  higher  pressure.   The 
assumption  is  made  that  slamming  will  still  occur  forward  of 
the  area  that  is  overestimated  and  that  the  pressure  can  be 
corrected  by  a  relative  angle  correction.   This  will  be 
discussed  later. 
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SLAMMING  PRESSURE  AND  HULL  FORM  COEFFICIENT  (k) 

To  estimate  the  slam  pressure  it  is  necessary  to  find  a 
functional  relationship  for  pressure  vs.  velocity.   As  mentioned 
earlier,  it  has  been  show i  that  equation  (2) 


k  r 


holds  very  well  and  can  bi  used  if  the  form  coefficient,  k, 
can  be  determined. 

Ochi  and  Motter  [16]  used  mapping  and  regression  analysis 
techniques  to  find  the  k  ralues  from  experimental  data.   This 
method  was  used  to  develo  )  a  computer  program  reported  in 
Reference  [16]  to  determi  le  the  coefficients.   This  program 
was  used  as  a  subroutine  ,n  this  work  but  it  did  not  give 
satisfactory  results  for  :>road  beam  sections.   The  regression 
analysis  had  been  performed  over  relatively  narrow  ship  sec- 
tions and  the  confidence  _ange  for  the  coefficients  was  0.073 
to  0.190  with  10%  risk  of  error.   In  applying  their  program 
two  problems  developed: 

1.  The  k  values  for  broad  beam  flat  sections  were 
well  outside  the  0.190  confidence  figure;  and 

2.  The  mapping  blew  up  in  many  of  the  broad  beam 
cases . 

Since  many  ships  of  comme  ~cial  interest  could  not  be  covered 
by  the  program,  another  m  ?thod  had  to  be  found  for  these  wide 
forms. 
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M.D.  Ochi  [17]  in  work  which  was  an  attempt  to  correlate 
the  k  values  from  two  and  three-dimensional  drop  tests  as  well 
as  model  tests  proposed  that  the  following  holds: 


B'2 
0.02  ^r  (6) 


where     B1  =  half-breadth  at  1/10  design  draft 

A'  =  half-section  area  at  the  1/10  design  draft 
This  work,  which  was  done  prior  to  the  regression  analysis 
technique  in  Reference  [16] ,  was  compared  to  the  results  ob- 
tained by  the  computer  program.   Better  agreement  inside  the 
confidence  range  was  obtained  if  a  slightly  higher  constant 
was  used.   Therefore  it  is  assumed  that 


B'2 
k  =  0.027  H— ,-  (6a) 


In  the  program  developed  in  this  thesis,  equation  (6a)  is 
used  when  section  values  above  0.190  are  obtained  from  the 
regression  subroutine.   This  is  done  since  equation  (6a)  yields 
larger  k  values  and  high  pressures  which  is  more  conservative. 

Now  that  the  form  coefficients  can  be  calculated,  it  is 
possible  to  determine  the  maximum  pressure  at  each  station 
along  the  ship  if  the  velocity  of  impact  can  be  found. 
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MOST  PROBABLE  EXTREME  PRESSURE 

The  impact  relative  velocity  is  approximated  by  the 
amplitude  of  relative  velocity  which  is  a  random  variable 
following  the  Rayleigh  probability  distribution  [14].   Using 
equation  (2)  and  the  Raylaigh  distribution  function,  the 
slamming  probability  distribution  function  is  given  by 


f  (p)  =  A  exp  [-A  (p  -  pQ)]  (7) 


p    -  p  <  oo 


where     p  =  slamming  pressure 

2 

p   =  threshold  pressure  =  kf . 

X   =  l/(k  R£) 

By  applying  order  st  itistics  to  the  probability  distri- 
bution (7),  the  distribution  of  extreme  value  of  pressure  in 
n  finite  number  of  impacts,  p  ,  is  given  by 


f  (p  )  =  [n  f  (p)  (F(p)}  n  1] 


=  nAexp  [~A(Pn  -  pQ)  ]  {l-exp[-A(pn  -  pq)  ]  }n  -1  (8) 


< 

p   -  p   <  °° 
^o    n 


where 

P 
F(p)  =  /   f(p)  dp 

Po 
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From  equation  (8)  it  is  possible  to  derive  the  most 
probable  extreme  pressure  in  n-observations,  p  .   Figure  3 
shows  the  probability  distribution  function  given  by  equation 
(8),   The  most  probable  extreme  pressure,  p  ,  can  be  obtained 
by  setting  the  first  derivative  with  respect  to  p   equal  to 
zero . 


d  f(P  ) 

a-p-2-   =  °  (9) 

n 


Using  equations  (8)  and  (9),  p   is  obtained 


.2 
'  * 


p^  =  k(r;  +  R^  ln(n) )  (10) 


This  most  probable  pressure  can  be  used  to  calculate  slam 
forces  and  is  one  of  two  pressure  options  in  the  computer 
program. 

EXTREME  PRESSURE  FOR  DESIGN 

Although  it  may  be  reasonable  to  compare  p   with  the 
largest  value  observed  in  experiments,  the  probability  that 
a  pressure  is  higher  than  p   is  0.632.   Since  the  probability 
distribution  function,  equation  (8)  as  plotted  in  Figure  3, 
is  never  zero  at  high  pressures,  the  extreme  pressure  for 
design  can  be  determined  only  if  a  small  probability  of  being 
exceeded  (a)  is  allowable.   The  extreme  pressure  for  design, 
p    (ot)  ,  can  be  found  from 
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i        f(pn»  d  pn  =  «  (11) 

pn 

Using  equations  (8)  and  (11) 


pn(a)  =  k  [±l   -  R£  :n{l-(l-a)1/n}]  (12) 


The  extreme  pressure  for  design  is  controlled  by  the  pre- 
assigned  small  probabili  -y,  a  ,  which  can  be  specified  by  the 
designer.   In  the  calculations  in  the  computer  program  Ochi ' s 
[14]  value  of  a  =  0.01  is  used,  but  can  be  changed  if  desired. 

SEA  STATE  INFORMATION 

In  the  calculations  necessary  to  determine  slam  pressures 
the  sea  state  must  be  specified  to  determine 

1.  the  relative  motions  and  velocities  and 

2.  the  duration  of  the  storm 

in  order  to  calculate  the  number  of  slams  from  equation  (4). 

In  order  to  make  the  application  of  the  program  as  simple 
as  possible,  the  only  pa  ameter  which  needs  to  be  specified  is 


the  significant  wave  hei'  ht,  H(l/3).   This  is  made  possible 
by  the  work  done  by  Louk  kis  and  Chryssostomidis  [15]  from 
which  the  r.m.s.  relativ  motions  and  r.m.s.  relative  velocities 
may  be  determined  for  th >  forward  stations  using  only  the 
parameters  of: 

1.  CD  =  Block  coef  icient 

2.  L/B  =  Length  to  Beam  ratio 
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3.   B/T  =  Beam  to  Draft  ratio 


4.  H(l/3)/L  =  S  =  Non-dimensional  sea  state 

5.  F  =  Froude  number  =  SHIP  SPEED/ /qh 

These  parameters  are  sufficient  to  give  good  initial  design 
values  for  relative  motion  and  velocity.   There  is  also  a 
dependence  on  longitudinal  center  of  flotation  and  radius  of 
gyration  but  it  is  not  particularly  sensitive  to  these 
parameters  [15]  .   These  tables  were  developed  from  Series-60 
hull  forms  and,  therefore,  can  only  be  used  for  ships  which 
can  be  approximated  by  th 3se  forms.   This  allows  calculation 
of  slam  forces  for  these  ships  without  running  seakeeping 
programs.   The  option  is  ilso  available  to  specify  the  relative 
motions  and  velocities  if  desired. 

The  storm  duration  £lso  needs  to  be  known.   Following 
Ochi  and  Motter  [14],  wav 3  records  have  been  analyzed  in  the 
North  Atlantic.   A  plot  cf  the  envelope  of  storm  duration  vs. 
significant  wave  height  is  shown  in  Figure  4.   This  envelope 
is  used  to  determine  the  maximum  duration  to  be  expected  for 
a  ship  operating  at  the  significant  wave  height.   This  envelope 
is  the  default  option  in  the  program  and  is  considered  the  ex- 
treme storm  for  design.   If  desired,  a  storm  duration  may  be 
specified  in  the  input  dc-  ta  to  the  program 
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FIGURE    4  -     STORM    DURATION    vs.    SIGNIFICANT     WAVE    HEIGHT 
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SLAM  IMPACT  REGION 

Using  the  time  of  st  orm  from  the  previous  section  as  well 
as  the  tabulated  r.m.s.  r itions  and  velocities,  equation  (4) 
can  be  evaluated  at  each  station  considered  to  determine  the 
number  of  slams  during  ar  extreme  storm.   The  number  of  slams 
is  greatest  at  the  bow  ar 1  decreases  rapidly  as  the  midships 
crossection  is  approachec  .   The  sections  where  N(t)  -  1  is 
assumed  to  be  the  longiti dinal  extent  of  slamming. 

From  model  tests  anc  drop  tests  [17]  it  has  been  found 
that  the  slam  pressure  is  essentially  dissipated  at  a  depth 
of  1/10  of  the  design  dr£ ft,  T1 .   Therefore  the  pressure  needs 
to  be  integrated  only  ovc  r  the  forward  portion  of  the  hull 
where  N(t)  -    1  and  below  1/10  design  draft. 

SLAM  FORCE  AND  IMPULSE 

The  information  is  r  ow  available  to  determine  the  most 
probable  and  extreme  pre£  sure  for  design  p  and  $  (a)  at  the 
keel  for  each  station  as  well  as  the  slam  impact  area  for  a 
particular  ship,  speed,  end  sea  state.  To  integrate  temporally 
and  spatially,  more  information  is  needed  about  the  spatial 
variation  of  pressure  in  the  girth  direction  as  well  as  the 
time  history  of  the  slam. 

The  local  pressure  1  ime  history  as  used  by  Ochi  [14]  is 
shown  in  Figure  5.   The  c uration  of  the  slam  pressure,  t. ,  is 
given  by 


t,  =  0.1  /L/520  =  4.4  x  10"3  /h  (13) 
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FIGURE    5   -    SLAM    PRESSURE    TIME    HISTORY 
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which  is  an  average  for  a  520  FT  vessel  whether  it  is  in  the 
flat  bottom  region  or  the  bilge.   Froude ' s  scaling  law  is  used 
for  ships  of  different  leigths.   The  slam  pulse  is  assumed 
triangular  with  equal  ris »  and  fall  times.   The  duration  of 
slam  is  one  area  where  it  is  felt  additional  work  needs  to  be 
done.   In  Appendix  A  an  a  .ternate  equation  for  t.  is  considered 
and  presented  here  as  equ  ition  (14). 


t±   =    3.6  x  10  3  B"     (seconds)  (14) 


where     B'  =  Half-beam  a   1/10  design  draft  in  feet 
Until  further  work  can  be  done  to  experimentally  determine 
which  time  duration  is  be  ;t,  either  can  be  used.   The  computer 
program  uses  eouation  (13   but  it  is  felt  by  this  author  that 
from  the  analysis  of  drop  test  results  beam  should  be  used  as 
the  scaling  parameter  ins  .ead  of  length. 

The  work  done  by  Och   assumes  that  there  is  no  dependence 
on  angle  between  the  keel  and  wave  surface,  and  that  the  worst 
slam  case  is  when  the  sla  i  pressure  travels  the  full  extent 
forward  or  aft  at  the  slowest  traveling  velocity.   It  is  felt 
that  the  wave  slope  has  a  large  effect  in  reducing  the  slam 
forces  predicted  by  Ochi ' s  slam  model.   Therefore  a  new  model 
of  slam  is  proposed  which  accounts  for  the  dependence  on  wave 
slope. 

Numerous  experiments  have  shown  that  the  slamming  becomes 
maximum  in  regular  waves  'hen  the  ship  encounters  waves  of 
length  L  on  the  order  of  1.0  to  1.1  times  the  ship  length,  L, 
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[12],  [13].   It  is  assume  1  that  a  similar  condition  has  been 
approximated  in  irregular  seas  leading  to  the  worst  slam 
condition.   For  the  relat  .ve  velocities  to  be  maximum,  the 
ship  will  be  pitching  dow  i  in  the  vicinity  of  zero  pitch 
angle.   For  this  to  occur   the  slam  will  take  place  in  the 
trough  of  the  wave.   Figure  6  depicts  this  model  of  slamming. 
The  slam  pressure  will  travel  from  the  forward  perpendicular 
aft  and  from  the  aft  slam  station  forward  meeting  at  what  is 
called  here  the  slam  cent  Br.   Assuming  the  wave  encountered 
has  height  equal  to  the  s  Lgnif icant  wave  height  allows  calcula- 
tion of  the  slope  between  the  keel  and  the  wave  surface.   This 
angular  dependence  will  b9  used  to  modify  the  slam  pressures 
from  what  would  be  expect  3d  for  flat  impact. 

Chuang  [9]  performed  a  series  of  drop  tests  with  varying 
deadrise  angles.   It  is  proposed  here  that  the  reduction  in 
pressure  with  deadrise  ar jle  can  be  used  for  similar  angles 
in  the  longitudinal  direction.   Figure  7  is  taken  from 
Reference  [9]  and  shows  a  significant  reduction  in  pressure 
with  deadrise  angle.   The  region  between  0  and  3  degrees  is 
where  the  compressibility  of  the  air  comes  into  effect.   It 
is  assumed  this  is  not  as  pronounced  in  the  longitudinal 
direction  due  to  the  escape  of  air  in  the  transverse 
direction.   The  higher  ar gles  are  almost  linear  when  plotted 
on  login  paper  and  the  slope  remains  almost  constant  for  all 
impact  velocities.   From  the  linear  region  of  Figure  7  a 
correction  equation  for  pressure  was  found  as 
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6  (a,0)  =  p  (a)  lO"0*12  f0 


n-  '  '    -n''  (15) 


where    £  (a,9)  =  maximum  keel  pressure  with  angular  correction 

Pn(a)    =  pressure  predicted  by  equation  (12) 

9       =  angle  between  the  keel  and  wave  surface 
in  degrees 

If  the  wave  is  assumed  approximately  sinusoidal  the  wave 
slope  can  be  calculated  as 


6  =  935.  H(l/3)  X  /  L  (16) 

where     X  =  distance  away  from  the  slam  center  in  feet 

Using  equations  (15)  and  (16)  the  most  probable  or  the 
extreme  pressure  for  design  may  be  modified  to  allow  for  other 
than  zero  impact  angle. 

The  peak  pressures  do  not  occur  at  the  same  time  along 
the  length  of  the  ship.   From  full  scale  measurements  on  the 
Wolverine  State  and  model  tests  [14],  it  was  determined  that 
the  pressure  pulse  traveled  from  forward  to  aft  or  vice  versa 
at  speeds  from  260  to  520  fps  for  a  520  ft  ship.   Ochi  assumes 
that  the  slowest  traveling  velocity  is  the  worst  slam  condition 
and  that  the  results  can  be  Froude  scaled  for  different  length 
ships  to  get 


vt  =  260  /L/52  0   =  11.4  /L  (17) 
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where     v  =  traveling  velocity  in  feet  per  second 
This  assumption  was  not  felt  valid  since  the  traveling  velocity 
should  increase  as  the  rate  of  entrance  into  the  water  in- 
creased— that  is,  as  relative  velocity  increased  so  should  the 
traveling  velocity.   It  is  agreed  that  if  a  slam  could  occur 
at  the  highest  possible  relative  velocity  and  the  slowest 
traveling  velocity  it  would  be  the  worst  case,  but  it  is  felt 
that  these  two  conditions  cannot  be  met  at  the  same  time. 
The  traveling  velocity  is  calculated  using  the  slam 
model  in  Figure  6.   Using  this  model,  the  height  above  the 
water  surface  at  the  slam  center,  the  most  probable  or  extreme 
relative  velocities  of  the  slam  center,  and  the  length  of  the 
slam  area  are  known  or  can  be  found.   From  this  information 
the  traveling  velocity  is  given  by 


vfc  =  1  /2t  (17a) 

t     s'   t 


where    1   =  slam  length  from  FP  to  aft  extent  of  slam 

t   =  time  slam  takes  to  travel  from  aft  extent 
of  slam  to  slam  center 


and 


t.  =  H/VEL  (17b) 


t 


where     H  =  height  of  slam  center  above  the  wave  surface 
at  start  of  slam 

VEL  =  most  probable  or  extreme  relative  velocity 
of  the  slam  center 
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The  last  information  required  before  integration  can 
proceed  is  the  spatial  distribution  of  pressure  along  the 
girth.   The  most  probable  and  the  extreme  pressure  for  design 
as  calculated  by  equations  (10)  and  (12)  and  modified  with 
equation  (15)  for  wave  slope  are  only  valid  at  the  keel. 
Following  again  Ochi * s  work  [14],  [18],  the  maximum  slam 
pressure  can  be  assumed  constant  along  the  flat  bottom  of  the 
hull.   From  the  end  of  the  flat  the  pressure  drops  off  approxi- 
mately linearly  to  zero  at  the  1/10  design  draft,  T' .   Figure 
8a  shows  this  for  a  typical  section.   The  pressure  component 
of  interest  is  the  vertical  component  on  the  hull.   The  values 
in  Figure  8a  must  be  multiplied  by  COS3  where  6  is  defined  as 
the  angle  between  the  normal  and  the  vertical  at  the  point  in 
question.   Figure  8b  shows  the  new  vertical  component  of  the 
pressure  distribution  needed  to  find  the  slam  force. 

All  the  information  is  now  available  to  calculate  the 
slam  force  and  impulse.   The  procedure  is  as  follows: 

1.  Using  equations  (4)  and  (5)  along  with  the  r.m.s. 
relative  motions  and  velocities  from  tables  [15] , 
calculate  the  nu nber  of  slams  at  each  station. 
The  last  station  aft  with  N(t)  -  1  is  the  longi- 
tudinal extent  o:  slam. 

2.  Calculate  the  most  probable  extreme  pressure  from 
equation  (10)  or  the  extreme  pressure  for  design 
from  equation  (12)  at  each  station  forward  of  the 
last  slam  station  as  found  above. 
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3.  Assuming  the  certer  of  slam  is  midway  between  the 
FP  and  the  last  station  where  slamming  occurs, 
calculate  the  slopes  using  equation  (16)  and 

the  corrected  piessures  at  each  station  using 
equation  (15) . 

4.  Calculate  the  vertical  pressure  components  at 

each  station  for  approximately  10  waterlines  between 
the  baseline  and  T1  as  in  Figure  8.   The  spatial 
distribution  of  the  vertical  pressure  is  now  known 
over  the  entire  slam  region. 

5.  Assuming  a  triar jular  shape  as  in  Figure  5,  and 
a  pulse  duration  of  t,  as  predicted  by  either 
equation  (13)  or  (14),  add  temporally  at  each 
location.   The  start  of  slam  pressure  rise  (t.) 

at  each  spatial  point  is  predicted  by  equation  (17a) 
in  the  longitudinal  direction  and  by  the  relative 
velocity  in  the  girth  direction. 
Due  to  the  numerical  integration  procedure  necessary,  a 
computer  program  was  written  to  perform  the  above  calculations 
and  numerical  integrations.   The  program  and  the  details  of 
the  above  procedure  are  discussed  in  the  next  section. 

By  evaluating  the  pressures  at  different  time  increments, 
the  force  time  history  can  be  determined  for  each  section. 
Integrating  these  force  histories  over  time  will  give  the 
slam  impulse  at  each  station. 
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V.   COMPUTER  PROGRAM  AND  DISCUSSION  OF  RESULTS 

A  computer  program  was  written  to  implement  the  procedure 
discussed  in  Section  IV.   A  listing  of  the  program  is  included 
as  Appendix  B.   There  are  many  input  and  output  options  which 
are  discussed  in  Appendix  C  which  is  designed  to  be  a  user's 
manual.   A  short  example  program  illustrating  the  data  deck 
and  some  of  the  possible  outputs  is  included  as  Appendix  D. 
in  this  section  the  results  obtained  from  the  program  will  be 
discussed  using  the  Mariner  and  Fotini-L  as  example  ships. 

Throughout  this  thesis  stations  are  numbered  from  0  at 
the  forward  perpendicular  to  20  at  the  after  perpendicular. 
Since  slamming  occurs  in  the  forward  portion  of  the  ship  where 
the  lines  change  rapidly,  half-stations  were  introduced  num- 
bered from  0  to  40  similarly.   The  slam  forces  calculated  are 
considered  as  centered  at  the  half-stations  extending  equal 
distances  forward  and  aft  of  these  half-stations. 

Two  computer  programs  were  written  and  used.   The  assump- 
tions and  calculations  oi  Reference  [14]  were  duplicated  and 
used  to  compare  with  the  results  in  this  work.   The  only 
deviations  were  the  use  of  equation  (6a)  for  k  >  0.19  and  the 
use  of  seakeeping  information  obtained  from  the  Reference  [15] 
tables.   The  results  for  the  Mariner  were  identical  to  those 
obtained  by  Ochi  [14].   The  significant  aspect  of  this  is 
that  good  results  were  obtained  without  recourse  to  seakeeping 
computer  programs. 
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OPERATING  CONDITIONS 

The  Mariner  was  used  as  one  of  the  example  ships  and  the 
computer  program  was  run  inder  various  operating  conditions. 
The  four  conditions  selec  :ed  for  discussion  here  are: 


1. 
2. 
3. 

4. 


MARINER 
MARINER 
MARINER 
MARINER 


FULLY  LOADED 
FULLY  LOADED 
LIGHT  DRAFT 
LIGHT  DRAFT 


7.4  KNOTS 
15.5  KNOTS 
15.5  KNOTS 

7.4  KNOTS 


all  at  a  significant  wave  height  of  25  feet.   Table  la  lists 
the  ship  characteristics  cor  the  two  loading  conditions  con- 
sidered.  This  ship  and  conditions  were  chosen  since  they 
were  used  in  Reference  [14]. 


MARINER 
BLOCK  COEF.  (C_) 

a 

L/B 

B/T 
DRAFT  FWD. 
DRAFT  AFT. 


H(l/3) 

LBP 
BEAM  (B) 
DESIGN  DRAFT 


FULLY  LOADED 

0.613 

6.95 

2.82 
27  FEET 
27  " 
25  " 
528  " 
76  M 
29.75  " 


LIGHT  DRAFT 

0.583 

6.95 

3.80 
16.4  FEET 
23.6 
25 
528 
76 
29.75 


TABLE  la  -  MARINER  OPERA: ING  CONDITIONS  AND  DIMENSIONS 
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A  second  ship,  the  F5TINI-L,  was  analyzed  and  is  one  of 
the  three  ships  reported  m  in  Reference  [4].  The  principal 
characteristics  and  loadi  ig  conditions  are  given  in  Table  lb. 


FOTINI-L 

BLOCK  COEF 

L/B 

B/T 
DRAFT  FWD. 
DRAFT  AFT. 


H(l/3) 

LBP 
BEAM 
DESIGN  DRAFT 
SPEED  (KNOTS) 


FULLY  LOADED 

0.84 
7.55 
2.79 
38.0  FEET 
38.0   " 
20  to  40  FEET 
800    FEET 
106 

38.0   " 
9.59,  14.38 


BALLAST 

0.65 
7.55 
3.55 
28.31  FEET 
31.46   " 
20  to  40  FEET 
800   FEET 
106 
38.0  " 
9.59,  14.38 


TABLE  lb  -  FOTINI-L  OPERATING  CONDITIONS  AND  DIMENSIONS 

R.M.S.  RELATIVE  MOTIONS  AND  VELOCITIES 

The  first  modification  to  the  Ochi-Motter  approach  was 
to  apply  a  broadness  factor  correction  to  the  r.m.s.  relative 
motion  and  velocity.   If  the  ship  responses  are  assumed  to  be 
narrow  band  processes,  the  r.m.s.  values  of  the  responses  can 
be  obtained  directly  fron  the  tables.   These  tabular  results 
do  not  take  into  account  the  fact  that  seakeeping  processes 
are  not  ideally  narrow  be nd  [15].   More  realistic  values  can 
be  obtained  if  a  broadness  factor  correction  is  applied: 
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mo  =  mo  (1  "  e2/2)  (18) 


where    mQ  =  r.m.s.  value  of  response  from  tables 

m'  =  corrected  r.m.s.  value 
o 


e   =  broadness  factor 

and 

2 

2         m2 
e   =  1  -  -4—  (19) 

mm. 
o  4 

where    m,  =  the  k   moment  of  the  spectrum. 
A  value  of  e  =  0.59  is  used  in  the  calculations  [7],  [15]. 
This  correction  factor  reduces  the  relative  motions  and 
velocities  and  therefore  reduces  the  slam  forces  predicted. 

RESULTS  USING  OCHI-MOTTEI  METHOD 

The  results  for  the  four  operating  conditions  using  the 
Ochi-Motter  method  with  the  broadness  factor  correction  are 
given  in  Tables  2  to  5  and  Figures  9  to  12.   These  results  are 
for  each  of  the  half-stations  where  slamming  is  predicted. 
Table  6  shows  the  total  slam  impulse  and  the  amplitude  and 
duration  of  an  equivalent  slam  for  the  ship.   Table  6  results 
are  presented  for  ease  of  comparison  of  the  two  methods. 
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MARINER         TOTAL  SINE  FORCE  APPROX.     TOTAL  IMPULSE 
CONDITION        AMPLITUDE       DURATION       (KIP-SEC) 

(KIPS)  (SEC) 


1  985  0.43  269 

2  4490  0.62  1766 

3  4406  0.68  1901 

4  1220  0.48  372 

TABLE  6  -  TOTAL  SLAM  FORCE  AND  IMPULSE  (OCHI  METHOD) 
NEW  METHOD  RESULTS 


The  changes  made  to  the  Ochi  method  discussed  in  the 
previous  sections  and  incorporated  into  the  second  computer 
program  are : 

1.  an  angular  correction  applied  as  per  equation  (15), 

2.  a  different  model  of  slamming  as  shown  in  Figure  6, 
and 

3.  a  traveling  velocity  calculated  by  equation  (17a). 
The  results  for  the  same  four  operating  conditions  are  given 
in  Tables  7  to  10  and  Figures  13  to  16.   Table  11  gives  the 
total  values  of  slam  force  and  impulse  corresponding  to 
Table  6  above. 


56 


•    ft    ft  ■ 


I  »-  — 

<  X   <_> 


•  ft  ft 

z 

o 

£  -•  — 

«  I-  U 

_l  «  Ul 

en  a-  it. 


Q  y  x  - 

Z3  Z  O  </) 

t-  -<a  a. 

•-  u>  a.  — 

_1  Q.     it 

0.  u.  <   — 

£  O 


ft    ft     ft   ft     ft 


_)  3    I 

:/>  a  a. 
r  — 


•   •    ■  •   » 


3  ruj  w 

r  <ua 

—  _iir  •-• 

X  U3X 

<  U.    — 

r 


•  •    •  •   • 

U.  i. 
ul  < 
X  -1 

r  o> 

z  u. 

c 

•  ft    ft  ft    ft 


D 

• 

a 

t/> 

• 

a. 

IT 

• 

— 

U' 

• 

a.  *— 

a 

• 

O   l/) 

Ci 

* 

tv  a 

• 

a.  — 

• 

V- 

■ 

►- 

X 

• 

(/) 

U. 

• 

o 

t 

r 

• 

* 

* 
• 

•   ■ 

• 

> 

7 

• 

»-  — 

c 

• 

•~  z 

ft— 

• 

u  *-• 

1  »- 

-    *  O  C 

O     •    _l     I 


»  _l  — 

ft  UJ 

*  a 


0.  U.'  — 
OC  It 
U  l-    — 


ft     ft   ft     ft 


ft    ft    ft  ft    ft 


»••■••••«>»■•■■■■ 

®  •   m  •   <j  ftoo  •     G  •<*) 

v»  •  »-<  •  fti  •  n  •   r\j  ft  •* 

s>  •  S  »  ri  •&  •  <s  •  & 

*  ft    •  ••  ft*  ft*  ft* 
SI  ft<9  •  is  ftS  ft    n  •  K> 

•  ft     »••»••     •■■•••■•• 

n        •(»)        ftr>        »  -r        •  .»        •  m 

•  *•        ft*         ft*        «•        ■    • 

S  •  S  ft  <S  ft  >S  ft    Si  ft  K 

•  ft      ••••••     («•■*»••> 

ro        t-t         *oj         »  <s        »cf>        »x 

•  ru        ftci        •  in        *■«-        •-# 

•  ■    ••••>•••••••••• 

■c  ■  ir         »  oc         •  -c        ftoo         •  ai 

O  *{T  *--i  ftp*  ft   -ft-  ft  0"\ 

*-i        •  p>        •*        •*        »  m 

•  ••••>••••••••••* 

•  ••        •*        ••        ••        •    • 

<n        »m       »(\i        »(v»       «(vj       •  .< 

>C  '    r^  *  s  •  tfl  ft   -h  «*-i 

*  n        •  •*?        •  oc        •  oo        •  r*- 

»    •    •••    ••    •    ••••*■    »•• 

»*.  »    fw  ft    C>  *  •*  ft    C.  ft   t*) 

i-»        i   c»        •  in         •  tu        •  ■ 

•  •••••••    >  ■    ■    t   *  •    ••• 

-*-oc*sCX>ftr^oc     •  -»    P~    •    .0  nj    *  K  in 

f;i~   »  ^  s.   •  c  <   •  r.  ex.  •  r».  ft.    •  .»  m 
•*        ftfnfvi«-trgft«*5f»>ft-#cvjftf*i^ 

•  *    .    •    •    •>•    •*•••*    »  »   • 

-*-*■   *ro&   i  h  >c   fts.<\i«ocoo   »Nif) 
i  cvj  p»    *<  (vj  r^   •«  rvj  si>  *•  nj  <>  «j  •-•  u%  *i  •■*  tf^ 
SS.ftS.lSftSlSftSSftS.5.     •  S.  S. 

•  •••••••••••••*•• 

6vb.»s.e.ftce.    •  &&••»&•■»,& 

•  ••>•••>>>•>■•    ••• 

J  »    Ifi  •   J.  •  -»  •   d  ft  u. 

fu        «ir        •■*)        •  r»        »o.        »-< 

S  ftS  ft   6.  ftSi  ft   S.  •   — 

•  ft*  ft*  ft*  ft*  ft* 
Cftf         *tk>         ■  c*          •  *-*          •  c£         •  s. 

•  >••*    >»•    •••«•*    •»• 
-t        •  fti        •  n         •«         •  in         •   * 

•  ■••■••••••••ft    ••• 


-I 

U 

J 
< 


a 

01 
O) 


Q 

o 

I 
»- 

UJ 

LU 

Z 


a 

UJ 

z 
a. 
< 


UJ 

_i 

CD 
< 


57 


•    •    •  • 


CO 

UJ  X  1-  - 
X.  <  Or  l_l 
•—    _J  "9    UJ 

>-  01 1-  in 
i/)  — 


o 

r  -•  — 

■«  i-  u 

_l  «  UJ 

it  or  in 

3  — 
O 


e  ut  *  — 

3  -Z  O  Ul 

I-  •-Q.  Q. 

•-  to  a.  >- 

_i      a  -* 

Q.  U.  <  — 

r  o 


•a  _i  oo 
_1  3    • 

to  x  a. 
x.  •-• 


t    i    >  •    f 


x.  <  u  a 

»-  jx  •- 


UJ 

• 

a> 

ct 

• 

UJ 

3 

» 

Bt 

IT 

• 

a. 

CT; 

• 

— 

u. 

• 

a 

a 

» 

c 

I/. 

u 

• 

cr 

Q. 

« 

u. 

— 

« 

i- 

• 

I- 

X 

» 

10 

UJ 

• 

a 

« 

X. 

ft 

1 

• 
• 

• 

ft 

• 

> 

7 

» 

ft- 

— 

C 

• 

•-» 

X 

» 

u 

i  ► 

•   It  c 

1  C 

O 

• 

_i 

1 

1 

• 

u. 

i. 

• 

> 

c 

_j 

• 

i: 

Uj 

* 

_l 

■— 

a 

• 

Uj 

1    u 

[r  u.  - 

C  O  :» 

u  t_>  — 


•     •     •    §     • 


X  4 

ft- 
10 

•     •     •    » 


•  •••••■ft  ftftft.       ftft.ftftft.i,.,,,,,,, 

Q        •  ■c  •  rvi  •  rv        •  m        ■  <j\  »  n        #rs        #ftJ        ,,£ 

G.         •  -*  •  r>  •  *         •   vc>         •  •>.  •  *©         «    j>        •  r*>         •  — < 

is         •  ■>  •  <s  •  s>         •  n.         »  e.  •  <s         •  15         •  s         •  s 

•  •     •  •      •  •     •  •      •  ft*  ft)     ft  ft      •  ft>     ft  ft      • 

IS  •   "S  ft    S.  ft   (3  »    s  ft   Si  ft    S  ft    Q  ft    ("2  ft    S) 

oj        •  ru  »rvj  •  ru        »(\j        •  m  •  r»>        #d        »  ^»        •  _r 

•-•        b  w  »  —  •>«-«        •  .-.        •  **  •  *■•        •  **        •«-«        »  -* 

•  ft*  •    •  »    •         ■    •         •    •  »    •         »•        «    •        #• 
st         •  '  •  i%  ■  »s         •  r.          •  ^.  *  *.         •  s        •  s         »  is 

•  •    »••»••  >   *    •    •    •   <     •  •    *  •  ••«■•••■•■• 

•  ft*  ••  ••         ••         B  •  ••         «    •         •    •         »    • 

--•  I   C  ftlfl  •  (V  9   D  •  r.  »CC  ft^  ft    X  •    » 

IP  •  oc  •   cr  ft  X  ft    r  •  -t  •  rvi  •   it  •   r-  •  *•< 

•  <vi  •-»■  •  r»         •  -tr         •  m  •  ru         •  x        •  o         «  jn 

•  •  •  ft—  tnj  ,nj  trt  i«  ■ 

ft   ft    •••  •  •   ■  •  •   ■    ft  •   ft    •••  •••    •  •••»■■»• 

a.         •  r*>  ••<  »rVj         ••<         •  r.  »  oo         •  .#■        •  ao        •-» 

D        ft(vj  #  x  •  *s         •  —        •  j~i  •  tri        •  -t        *(*>        •  oj 

•  c\*  •   n  •   x         ••-*         »  cr>  »  x  •  X)         »  <\t         »  x» 

•  ft     ft    ft     ft  ft    •     ft  >•••••     ft   ft    •  ft    ft     ft     riftftftftftftft 

•  ft*  •*  ft*              ft      •              ft     •  •     •  9ft  ft       *  ft      • 
ftv             ft.30  »(\J  •    -t             »X             »-H  »rv             #r»            •■»             •►» 

t^        •(»)  «j»  •  jj        <<h        ft  j  #ijn        »m       ■  <\j       »-< 

•  ■»  «r^  »^        »»«        «r^  *ir        »o>       •<*)        »m 

•  •  #*-•         »<\j         ftfn  »r)         »c\j        »M        ••-• 

•  •    ■••  •   •    •  »  »    •    •••    »»•  •••    ■••■*»•»■ 

•  ft*  •    •  ft*         ••         #•  ft*         •    •         ■    »         •    • 
s,         •-•  »r-  ftc>         »(v.         •  u.  •  u         »if.        »«r         »(\i 

IT.  ft    —  •&  •-#  •    —  ft   C\j  ft    Jj  ft   <M  •  • 

(T        «i~  »if>  »f-        «r\.        »*-■  •  •  •  • 

•  •■*••••  Sftftftftftftftllftftftftftftftftftftl 

r^-rvj   i  ac  <£  ftioa  •*>©   *  if  n   »n  o^  •  aj  -r    ftD^^ajrv   »s:<r 

•  •ft**ft**ft*****ft**ft**ft**ft**Sftft 

<  t.   «  >■  f  •  f  f  •CN<a«>ii'4  •ocoL»oc(T)»otir»-»>-« 

*i         •  -•    (\.  •  ur.  roftot*    #-.*•*  a,  ftCvj.t,ftoc-*ftir;(Vjftn*-« 

•  ftftftftftftft  ftftftftftftftftft  ««*     ftftftftftft     ft   ft    ft 
a*  a.  ftxui  ftr*-*  •  \s.   l    ft^  ro  «cwuc  »*^r*  »c^c  •  oc  tr  •  -o  s 

(Vft  0s    •<  fVJ  O^   *<  (Vft  0^   «•  fu-  0*J-   *«  ra.  CO    •«  (\i  h*     ft*  «MI  tx-   *1  *-l  *»}    ••  *H>  •*>   o  -^  >0 

&  1    ft  &  ii  •    ...    ,..  ftise.    •   r,    .•>    •  i.    ,-s.  •  n    c>    ft&is    ft  c>   (S.    ft(S>S 

•  ••••ft**ft*      •*•••••*••#••■••••• 

cj>c>,ftuts  ftCt>c.  •  ^  &■•&.&.   •  c*  c.  *4SCs.ft6kC.ftCC>*c.e. 

■  ft    ft  ft    ft  ft   ft    ft  ft   ft    ft    ftftftftftft  •■•»■*•*•■•• 

<ft  •    a  •  If  J  •    -»  •    r.  *u  *    0  •   l^  ft    -jL  «<jjj 

o.  ft    u  •    s6  ft>r^  •    cr  •    —  •    r  •    ir  •    r^  «    x 

&  ftC»  ftS.  •&.  »5».  •—  •-■  «—  ••-  **-• 

•  ••  «•  ft*           ••           ••  ft*  ft*  •     •  •      • 
C            »Sj  ftC  »6.            •&.            ftS  *6>            •&,           ftS.  »&>> 

■  ft      ••■  •••  ■•«      •«•      •••  >*>«■••■•••• 

*<        »c\)  •(•>  •-»        »ir        ft-o  •*»        »oo       •  •>        •>*< 

•  •    •••  •■•  ••■    •■•    *••  •••    .••••.    »«• 


3 
U 


cr 
3 
VI 
If) 


a 
o 

x 
»- 

UJ 

Z 
$ 

UJ 


or 
UJ 


a. 

< 

2 

i 

CM 


GO 

UJ 

_l 
CD 
< 


58 


■    ■    •   •    • 


HI    1     I-    - 

r  <  a  o 

-    J«UI 


u- 


—  w  a.  — 
_l       xv 

B,  la.  <   — 

t  o 


111 

<_) 

i:  o- 

UJ 

<  _i 

i/j 

_l  3 

• 

cn  a 

a 

r 

— * 

•-• 

■x 

ft     ft      ft    •      • 


3  r  uj  in 

i  <  (_>  a 

>•  ja  » 

x  to  cj  * 

<  u.  — 

J- 


a  jl 
uj  < 

O)  _l 

r  in 

3 
Z  U 

o 


_     •  - 

a.    •  u. 

IT     •  0 
(/■'ft  — 


q    •  q   a 


•  > 

Z      »    t- 


c    • 


>-  1 

»  l_(  — • 

'  ►-  •  r  r  i  r 

C     ft    ^  I 

1     •  u.  z 

•   >  c 

-J   •  z 

U      •    _J  — 

J     •    . 


•     •     •    •     • 


I  u. 

a:  uj  — 

c  c  v 

„.   O  — 


•     •     ■    ft     • 


I  «1 

ft- 


•  •    ft  ft    ft  ••    •  I    •    •    >    .    ■  •    •    ■  .>.»........,,, 

Si  ft   CC  ft    J)  ft    C)  •   n  #(7>  »^,  ,    ™,  #^,  #c0  #s 

c>         •  — •  »  ri  »ir  »rv  •  oc  •  rv        •  ;n         «  n        »  „        #s 

S         •  S  ft  <S  ft  s  •  s  •  s>  •  s        •  s>         «<s        >  s         •  <s 

•  »    ft  ••  ••  ••  ••  ••         ••         »    .         •     .         •    . 

S.  •   S  •'S  ft    <S  ft   S.  ft  <S  ft  K  ft  S  ft    S  ft    S  ft  O 

••••••»»•»••••••»»•»••••»»•»»»»» 

T*  «-H  ft    **  ft     ^  ft     ^  9    ^  t    ^  #^,  #^,  ft    ,«•  ft    ... 

(\i        »nj  »<\j  •  cu  »(\j  •  r>  »n        •  r>        •  .#•        •  ■*        •  in 

•"•  •     ■-•  ft     — «  ft     «-«  •      r-.  ft      .-  ft    r-4  •      ^  I      rt  ft      «*  ft     V* 

•  •      •  •     •  •      *  ft      •  ft      •  ft*  ft      •  V*  ■•  ft     • 

fi        •  s  »  is.  •  is  •  s  »  r  •  is        •  c         •  s        •  rs        •  s> 

•  ••••■•»■•■••**•••»•*•«••■•■•■»• 

r*         •  cc  •  -*  •  *c  •   .*■  •  io  *>£         •  r*  »  tr         •   ao         •(VI 

-»•        •  ~o  •*  »(\j  •  .»  «s  *  n        •  i^        »  .»       •  oj       «ip 

•  tv,  •-»  in  •  n  »{\j  •  —        »r^        #p)       *(n        «4) 

•  •  •  •.*»CVJ*(Vj         ••«••«•  » 

•  •    »  •   •  •  »   •  ••>  »   •    •  •  •   •  »•••»•    ••■•••••• 

-o        •*  »m  »(r  •  j-  #o>  •«       »(vi        •  rt       •♦       ••«. 

•  **  V*  V*  »      •  ff     •  •     *  •      •  ft*  ft*  •     * 

o        »s  ftir.  »uo  ftm  »m  •*        •«        ••<       •<•)        »fg 

•  (Vj  •  m  sun  *s  »  a.  •  r-         •  ■*         fttvi        •  «        •<« 

•  •    •••  »  »    •  •••  ••»  •••  •  •    >    ft   ■    •    *•••••»** 

*         ft*  ft*  ft*  ft*  ft*  ■•         ft*         ft*        ft*        •• 

r\j        ft'ji  •  ^  »fv  »s  •  so  •*»        »x        •-»        •-»•        »r» 

f*        »s  •  <r\  ftin  •*«.  ••-•  •  vo        •*»>        »r»        •*        ftio 

•  -r  <<C  *S  •  o-'  ftuo  •  m         •  ac         ftCVJ        fttf         •   " 

•  •  ft^i  ft—  »ri  »rj       ftcvj        »w       ••<       •*< 

•  •  ■»•  >  »  •  •■■  *••  •■•  •••••■  •••••»  ••• 

•  ft*  ■•  ••  ft*  ft*  ■•  ft*  ft*  ft*  ft* 
<r  •&  ftD  •  "u  •  r,  •  a  •  a^  **^  to  •  -»  •  oj 
•0         •  c.  »<\j  •    -c  •   ^t  «r)  ttCi         •   -r  *  IT'         •«-•         • 

r^         «(V  »rv  »cvj  »x  ftuo  »rw         ••-<         •  •  • 

<\Jft(Vjft.-ift<-ift  •  •  •  •  •  • 

•  ft    •■•  •••  ••<  •••  ••>  *»••»    ■    •»••»•••• 

nk    t  n  p  ftirrv  •  -o  ?s  •  r\j  is  »-r-»  ftf^wftpvs   t^st-ftiDft^<o 

•  *ft**ft**ft**ft**ft**ft**ft**ft**ft**ft** 
f^c    *  a   j  t  /   ^  i   j    «  t  ■,    t.  ft  a   a  •  --  -?    «    -»    r .    »^rv«-*r^ft<v<j 

—  •fT>cvjftiT-rri  ftr^*  ft^-.r^  •  a"    a  »fvjr^*otuo    i  tfi  n  •  nn  ifti 

•  •     •••  •»■  ■    •    «  *  #    «  «■•  «  »    «     ftftft«««ftft«ftftft 

■0-*»uorv  ».*(**  ft<vx**-i-»  »a<G  ftoouo»^-^ft^^ftri(\*«ft.ot 

i\j  ctn    •<  f\j  oo-  ••  nj  oO-  •*  f\>  r*    himn  urtr    •>  «-i  »0  fti  ^^  ^  •!  «r*  m  •!  •—  trv  •>  *-»  -• 

'-    ' .     •&%  ft  ?i    &  ■&&-  ft   C£   &  •    '-   <~--  ■   S  G     •    ^-    '■     «GGftC»S«SG 

•  *ft**ft**ft**ft**ft**ft**ft**ft**ft**ft** 
C*    &     ftCkCk  ft    £*    r_.  I    c    C*  I    a   &  I    C    b  I   L   &    •    &   S     «6&ICGftS& 

•  ft      ••••••■■•••••>•■••      ■■•■•*••■•■• 

-»        «iri  «ir  ••»  »n  «in  •»>       •  i"-        •*       «-c       •■* 

(V.  fttl  ft**  »r^  •    <f-  ••-<  •(*>  ftlf>  ft*"^  ftflt  ft(VI 

G  »i  •&  •&  •—  »--  ••-  *  •*  ■•H  •-*  •   >M 

•  *•  ••  ft*  ••  ft*  ft*  ft*  ft*  ft*  ft* 

e>        •  -^  •  r,  »  r,  «g  •  "*  •  r-        •  r-        •  s       is       ftSi 

«*    •»•  •••  <••  ••»  »••  •»••••    ••••»•••» 

—  »(Vj  »n  •-»  •*/)  ft^o  »r^        ftac        *  9\       »ts        •«-• 

*■  ••■  •>■  •••  •••  •••  ••••••  ••••••••• 


to 

z 
o 


< 
-I 

3 


a 
o 

i 
i- 

UJ 

Z 

Z 


cr 

UJ 

z 
or 
< 

i 

CO 


<n 

UJ 

_i 
CD 
< 


59 


in  ri-  - 

I   <i  u 

"    J4  U 

►-    {/)  ►-  05 

IT  — 


r  —  — 

«  t-  u 

_>  <  u. 

u  a-  o 


UJ 

O  UJ    X  — 

t-  «a  a 

•-.  en  3.  ►-. 

-i      a.  * 

a.  u.  <  — 

t  o 


<  _i  i/i 

_i  3    I 

i/>  a.  a. 

jt  ■-• 

-<  x 


3  [y  i/i 

2.    <i  u  a 

<       u.   — 


</j 

Uj  < 


IT    •  0. 
1/-,    •  — 

UJ     ft  IX    *-i 

a    •  o  (jo 

a    •  tt  a 

»  x  — 

■ 

►-   »  ►- 

x     •  (/> 

u.   •  o 

•  £ 


•  >- 

»   ►—    — 

»  ■-  r 

•  o  —• 
•   C  '  c 


•j-  u 

a  uj  — 
u  u  — 


•  •     »    ■     • 

u  C 

_l  •— 
<  *- 

I  < 
in 

•  »    »   •    » 


-••»••••••••••    •»•••• 

is         •  in         •  cr\         •*         •  tr>         •  ul  «s 

~«        »  —        •  w        •  -t        •  <\j        ■  .-.  »  c 

s  •  r.  *s  •  Gj  >  s  •   rj  »  q 

»  ft    ■  »    •  ft    •         ft    •         •*  »• 

<S  •    Q  »(S  »«  ■   s  ft®  ftlS 

»»■■••«•••••■■■■••■• 

•"•  •     »-•  *    »-«  ft*-*  ftv-4  ft     — •  ■  *-• 

m         •  c>         «r)         »  m         •  -»         •  .#  •  j) 

»  ••  ft*  ft*  «•  ••  ft* 
a         •  «        •  s.        »  e»        •  s.        •  s.         »c« 

•  ft*  ft  •  *•  ft*  ••  #• 
-»  *3  »  <"  fttv,  •  .jo  ftm  »  r» 
-»           •-«•           •   a>           ft  0-           »(Vj           •   r^  .  rv 

•  <vi         »  en        •  ir>        »m         .  _r  fto 

•  ft    •••    •»»••••■•    •••  ■•» 

•jo        •  s        »oo        »p>        •  o        »ir>  •  ro 

•  ft  *  *•  ■•  •>•  ft*  ft* 
m        »o        »fo        ftcr         •  r»        »<vj  »  *o 

•  o.  •  ?>  ft*  ft*  ft  •#  ft  P> 

•  ••»••»•••••■•>•>••• 

-<  ft*  ft  G  •  •*  »    sC  ft   <\J  »». 

i-»        »oc        »(\i        ftM        •  ir>       •►»  »ir> 

•  <*)  ft>0  ft  0»  ftCC  ft  Is*  ft  >*) 

ft*    •••    ftftftftftftftftft    *••  ft  ft    ft 

*-        «w         »c         •  r>.        »  re        •  f»  »* 

-»■        •—        »   -c        •  <r        •  x        •  (v  • 

a        •  -o        •  m        » -i        •  ■  • 

■  •    ft  ft   #    »»••••«»»    •»•  ft  »    * 

r*cc    »xrv    ft(\i(TftoenjftmeC'ftcsao  •  is  *-■ 

•  *ft*  •■*  *ft**ft**ft**  ft*  • 
ex  isri-r'iiftKia  *  >c  G  •  ir  <»> 
»-•        •nfv.»»ofT>#'*)-»«incvjftrr)(\jftfti*-< 

H«ft*ft«ftftftlftftftft*llfl  ft    «     « 

i\,  c**  •  *-•  ir.    »?<r4ixNt<cn   •  in  x  •  (*}■* 

i.   E.    •    G.   K    •  G    S    •   G    £.    •  C   5.     »ss  »cs 

•  •   ■•••••■••■••■••••* 

■  ft      ftftft      ■■■«■•■■■      •■■  ft    ft     ft 

•  •  ir  «ir  «  *  »r».  »  ur»  •  <r 
n»  »it  •  *c  •  r»  •  0"  •—  »«*) 
g        •  s        »<s        •  x        •  <s        •  —  »•-• 

Si  »    -..  »  &  •    5.  •   r,  •    ■:  >S 

•  ■  •■•  »»••»••••  •>•  ftftft 
-~            »<\j            •D            •-»             •!/>            •**)  »^ 

•  ■•••■•■•••••••■•••ft 


o 


cr 
3 
IT 


Q 
O 

I 
f- 

uu 

UJ 

z 


a: 

UJ 

z 
< 


UJ 

_i 
CD 
< 


60 


HI 

6 


HI 

M 

•- 


m 
6 


■ 
o 


o 


Ll) 
Q 
O 

LlI 

z 


q: 

LU 

Z 

< 

2 


i 

LU 

tr 


(SdlW    TECH 


61 


in 


m 
o 


I 

a 


fU 
d 


ni 


8   §   § 


o 

d 


UJ 
O 
O 

LU 


UJ 

z 
< 


UJ 

cr. 


(SdIM)    TWJ1 


62 


UJ 
Q 

o 

UJ 


(Z 
UJ 

z 
a 
< 

i 

CO 


UJ 

q 

u 


(SdIV) 


63 


i 1 1 1 1 r 


i r 


in 
6 


m 

o 


8 


J 1 1 1 1 L 


J I 


UJ 
Q 
O 

z 

UJ 


UJ 

z 

< 

2 


UJ 

or 

U. 


(SdXW 


64 


MARINER  TOTAL  SINE  FORCE  APPROX.     TOTAL  IMPULSE 

CONDITION        AMPLITUDE       DURATION         (KIP-SEC) 

(KIPS)  (SEC) 


1  1632  0.18  187 

2  7312  0.21  978 

3  7105  0.22  995 

4  2084  0.17  232 

TABLE  11  -  TOTAL  SLA 4  FORCE  AND  IMPULSE  (NEW  MODEL) 

COMPARISON  OF  RESULTS  BETWEEN  THE  TWO  METHODS 

In  general  the  new  model  predicts  lower  slam  impulses, 
shorter  slam  durations  an  1  higher  slam  forces.   In  addition, 
the  time  histories  are  mo  "e  triangular  in  the  new  model. 
The  lower  slam  impul 3es  are  a  result  of  the  angular 
correction  term  only.   If  the  correction  were  not  used,  both 
methods  would  yield  the  same  impulse  values. 

The  shorter  total  duration  is  a  result  of  the  slam  model 
and  the  calculated  traveling  velocity.   The  model  reduces  the 
duration  by  approximately  1/2  since  the  pressure  wave  travels 
forward  and  aft  toward  the  slam  center.   An  additional  decrease 
in  slam  duration  occurs  because  the  calculated  traveling 
velocities  are,  in  general,  higher  than  those  calculated  by 
Ochi.   For  these  four  operating  conditions  the  traveling 
velocities  range  from  731  fps  for  case  3  to  1016  fps  for 
case  1. 
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The  total  force  on  the  ship  is  the  spatial  and  temporal 
integration  of  the  local  pressures.   In  the  new  method  the 
pressures  are  always  less  than  or  equal  to  the  Ochi  method 
but  the  new  method  yields  larger  forces.   This  is  due  to  the 
integration  over  time.   Referring  to  the  plots  of  both  methods, 
Figures  10  and  14,  the  reason  for  this  can  be  seen.   The  half- 
station  forces  are  more  bunched  together  in  the  new  method. 
The  total  force  on  the  ship  at  any  time  is  the  sum  of  the 
forces  on  each  half-station  at  the  same  time.   Due  to  the 
greater  overlap  in  the  new  method,  more  half -stations  are 
summed  at  the  same  time,  resulting  in  higher  forces. 

The  importance  of  the  various  parameters  of  storm  duration, 
threshold  velocity,  local  slam  duration,  significant  wave 
height,  and  most  probable  vs.  extreme  pressure  for  design 
were  looked  at  to  determine  the  sensitivity  of  the  results 
to  their  variation. 

SENSITIVITY  TO  STORM  DURATION 

To  determine  how  important  accurate  prediction  of  the 
maximum  storm  duration  was,  the  MARINER  case  2  was  run  for 
different  values  of  storn  duration  with  all  other  input 
variables  held  constant.   As  seen  in  Figure  17  the  total  slam 
impulse  is  not  sensitive  to  storm  duration.   A  50%  increase 
in  storm  duration  from  3C  to  45  hours  results  only  in  a  4% 
increase  in  slam  impulse. 
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SENSITIVITY  TO  THRESHOLD  VELOCITY 

Since  the  threshold  velocity  at  which  slamming  occurs 
has  been  defined  differently  by  various  authors,  the  sensi- 
tivity of  slam  impulse  to  changes  in  threshold  velocity  was 
determined.   Values  of  threshold  velocity  between  6  and  20 
ft/sec  were  used  for  MARINER  case  2  and  the  total  slam  impulse 
plotted  against  threshold  velocity  in  Figure  18.   From  this 
figure  it  can  be  seen  that  the  use  of  12  ft/sec  is  reasonable 
and,  at  worst,  conservative. 

SENSITIVITY  TO  LOCAL  SLAM  DURATION 

The  total  slam  impulse  on  the  ship  is  the  sum  of  all  the 
local  impulses.   Those  parameters  which  affect  local  impulse 
also  influence  the  total  impulse  the  same  amount.   From 
Figure  5  the  local  slam  impulse  is  given  by: 


h   -   l'2   Pn(a)  fcl  (20) 


2  2 

where     Ip  =  local  impu] se/in   =  lb-sec/in 

From  (20)  it  can  be  seen  that  a  large  variation  in  t.. 

results  in  a  large  change  in  I.   and  subsequently  in  the  total 

impulse.   Therefore  total  slam  impulse  is  very  sensitive  to  the 

local  slam  duration.   Sirce  the  determination  of  t1    is  in 

question,  this  is  considered  the  area  where  further  work  is 

necessary  (Appendix  A) . 
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AFFECT  OF  SIGNIFICANT  WAVE  HEIGHT 

Another  ship,  the  Fotini-L,  was  also  analyzed.   Several 
computer  runs  were  made  using  both  programs  varying  the  signi- 
ficant wave  height  from  12  to  40  feet.   Figure  19  is  a  plot  of 
the  results.   The  Ochi  method  gives  results  approximately  2  1/2 
times  the  new  model  results.   As  expected,  the  impulse  in- 
creases with  significant  wave  height  and  with  velocity  in  the 
seaway.   The  rate  of  increase  decreases  at  higher  significant 
wave  heights  and  is  expected  since  the  ship  responses  behave 
similarly  at  higher  sea  states. 

MOST  PROBABLE  vs.  EXTREME  PRESSURE  FOR  DESIGN 

The  four  Mariner  cases  were  run  using  the  new  method  for 
the  most  probable  extreme  pressure  and  the  extreme  pressure 
for  design  cases.  Tables  11  and  12  summarize  the  results. 


MARINER  TOTAL  SINE  FORCE  APPROX .       TOTAL  IMPULSE 

CONDITION  AMPLITUDE      DURATION  (KIP-SEC) 


1  733          0.20  91 

2  3248          0.25  510 

3  1-663          0.25  591 

4  1147          0.19  137 

TABLE  12  -  TOT/L  SLAM  FOi CE  AND  IMPULSE  (MOST  PROBABLE 
EXTREME  PRESSURE  USED) 
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The  duration  of  the  slam  is  greater  in  the  most  probable  case 
since  a  lower  relative  velocity  predicted  results  in  a  slower 
traveling  velocity  (n.b.,  Table  11).  The  impulses  are  approxi- 
mately 1/2  those  obtained  from  the  extreme  pressure  for  design 
calculations.  Since  the  most  probable  extreme  value  is  indi- 
cative of  the  largest  value  observed  in  a  trial,  these  values 
should  be  used  if  any  comparison  is  done  with  full  scale  data. 

SHAPE  OF  THE  SLAM  FORCING  FUNCTION 

Knowledge  of  the  shape  of  the  slam  forcing  function  is 
necessary  to  input  the  slam  force  into  the  S.H.V.R.S.  program. 
The  S.H.V.R.S.  program  will  accept  three  input  shapes  for  the 
forcing  function: 

1.  triangular 

2.  half -sine  wave 

3.  damped  half-sine  wave. 

The  forcing  function  can  be  applied  at  each  of  the  half-stations 
with  time  delays  corresponding  to  the  traveling  velocity  or  a 
single  slam  force  can  be  applied  at  one  location.   Since  much 
of  the  previous  work  [2] ,  [3] ,  [4]  used  a  single  sine  pulse, 
this  program  develops  a  half-sine  approximation  to  compare  with 
these  results. 

The  shape  of  the  force  at  each  half-station  is  different 
in  the  Ochi  method  and  this  method.   In  the  Ochi  method  the 
assumed  shape  of  the  local  slam  (Figure  5)  is  not  critical. 
Several  pulses  of  different  shapes  were  tried  in  the  Ochi 
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program  but  the  predicted  half-station  forcing  functions  were 
about  the  same.   This  is  because  the  traveling  velocity  is 
slow  enough  that  the  slam  starts  at  different  times  along  the 
length  of  the  half -station.   In  the  new  method  the  traveling 
velocity  is,  in  general,  high  enough  that  the  slam  starts  at 
the  same  time  along  the  entire  half-section.   This  results  in 
a  force  shape  the  same  as  the  assumed  local  slam  pulse.   The 
shape  of  the  slam  force  at  each  half -station  (Figures  9  to  12) 
in  the  Ochi  method  can  be  approximated  as  half -sine  functions. 
The  program  calculates  the  amplitude  and  duration  of  these 
approximations.   The  shape  of  the  force  at  each  half-station 
(Figures  13  to  16)  in  the  new  method  can  be  assumed  triangular, 
of  height  equal  to  the  maximum  force,  and  of  duration 


t,  =  1.1  t,  =  4.8  x  10~3  /h  (21) 


where  t„  was  determined  by  neglecting  the  "tails"  of  the 
force  function  and  maintaining  the  same  impulse. 

The  total  slam  force  shape  can  be  obtained  by  adding  the 
forces  from  each  half-station  at  the  same  time.   Figure  20  is 
a  plot  of  the  total  slam  force  for  the  MARINER  case  2,  time  of 
storm  30  hours,  calculated  by  the  new  program.   Figure  21  is 
the  plot  of  the  forcing  function  at  each  half-station.   If  the 
points  which  are  less  than  1%  of  the  maximum  value  are  neglected 
and  a  linear  regression  analysis  done  on  the  remaining  points 
in  Figure  20,  the  triangle  formed  is  almost  equilateral.   There- 
fore this  shape  should  be  used  for  the  total  slam  forcing  function 
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MIDSHIP  SLAM  BENDING  MOMENT  AND  STRESS 

The  FOTINI-L  was  considered  at  different  significant 
wave  heights,  velocities,  and  loading  conditions.   The  results 
for  the  total  slam  impulse  in  the  ballast  condition  are  shown 
in  Figure  19.   As  the  significant  wave  height  and/or  ship 
speed  is  increased,  the  slam  center  moves  aft.   Donovan  [4] 
determined  midship  bending  moment  and  stress  for  the  FOTINI-L 
using  a  100  ton-sec  half-sine  slam  impulse,  applied  36  feet 
aft  of  the  forward  perpendicular.   This  slam  location  is  close 
to  half-station  2  of  this  model.   From  the  computer  runs  on  the 
FOTINI-L,  two  operating  conditions  resulted  in  slam  centers  at 
half-station  2.   These  are  given  in  Table  13. 

CONDITION   LOADING     SPEED     H(l/3)    TOTAL        DURATION 
CONDITION   (KNOTS)     (FEET)   IMPULSE      OF  IMPULSE 

(KIP-SEC)        (SEC) 

1-FOTINI-L    FULL      14.38       32      1199  0.18 

2-FOTINI-L   BALLAST     9.59       20       960  0.18 

TABLE  13  -  FOTINI-L  -  SLAMS  PREDICTED  AT  HALF-STATION  2 

From  the  graphs  of  response  in  Donovan's  work,  these 
slam  impulses  give  the  moments  and  stresses  in  Table  14  for 
various  stiffnesses.   This  is  determined  assuming  a  linear 
relationship  between  bending  moment  and  impulse  which  holds 
if  the  slam  durations  are  the  same — that  is,  holding  all  else 
constant,  doubling  the  slam  impulse  doubles  the  bending 
moment  response. 
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BENDING 
MOMENT 
(FT-TONS) 


CONDITION  60% 
1-FOTINI-L  343 
2-FOTINI-L    246 


STIFFNESS* 
80%     100%     120' 


364 
268 


375 
283 


391 
296 


140% 

404 

313 


STRESS 
(KPSI) 


1-FOTINI-L    26.5    22.1    18.9     16.5     14.5 
2-FOTINI-L    21.9    19.2    16.7     14.5     13.1 


TABLE  14  -  FOTINI-L  -  BENDING  MOMENT  AND  STRESS  AMIDSHIPS  DUE 
TO  SLAMS  CENTERED  AT  HALF-STATION  2 

*The  stiffnesses  are  %  of  as  built  values. 


The  values  in  Table  14  for  the  fully  loaded  case  are 
approximately  5.4  times  what  was  predicted  in  Reference  [4] 
and  4.3  times  the  predicted  results  in  the  ballast  condition. 
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VI.   RECOMMENDATIONS  FOR  FUTURE  WORK 

As  seen  by  the  values  of  stress  obtained  in  Table  14, 
the  slam  forces  and  impulses  result  in  large  midship  stresses 
on  the  FOTINI-L.   The  stresses  could  be  evaluated  only  for 
slams  at  half-station  2  since  dynamic  response  results  were 
available  only  at  this  station.   As  the  slam  center  moves  aft, 
the  total  slam  impulse  and  force  increases.   Higher  slam 
stresses  could  result  from  slams  at  different  locations. 
Therefore  additional  slam  locations  should  be  considered  in 
the  S.H.V.R.S.  program  in  any  study  of  the  slam  response. 

The  shape  of  the  total  slam  force  is  more  triangular 
than  sine  (n.b.,  Figure  20).   The  S.H.V.R.S.  program  should 
be  run  for  several  triangular  pulses  and  the  results  compared 
to  the  half-sine  results  available.   If  they  do  not  come  close, 
a  triangular  pulse  should  be  used. 

The  validity  of  using  a  single  total  force  approximation 
vice  distributed  forces  along  the  bottom  needs  to  be  verified. 
The  results  of  this  work  give  the  distributed  forces  as  well 
as  a  half -sine  approximation  for  the  total  slam  force.   Both 
should  be  run  on  the  S.H.V.R.S.  program  to  determine  if  the 
two  agree.   If  they  do  not,  either  the  distributed  forces 
should  be  used  or  a  new  approximation  developed  which  does 
agree. 

Although  the  slam  impulses  are  less  in  this  method  than 
in  the  Ochi  method,  they  still  result  in  large  slam  stresses. 
There  is  still  reason  to  believe,  because  of  these  large 
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values,  that  these  impulse  values  are  also  overpredicted.   As 
discussed  earlier,  and  in  Appendix  A,  equation  (13)  for  the 
duration  of  the  local  slam  pulse  is  in  question.   There  pre- 
sently exists  a  large  amount  of  literature  on  drop  tests  but 
few  model  and  full  scale  results  are  available  to  determine 
durations.   Correlation  with  full  scale  data  needs  to  be 
accomplished  to  determine  if  equation  (13)  is  valid  or  some- 
thing like  equation  (14)  would  be  better. 

A  series  of  model  tests  should  be  run  to  verify  that  the 
model  of  slamming  used  here  is  valid.   This  would  include 
verifying  or  modifying  the  traveling  velocity  calculations 
and  the  correction  equation  for  the  angle  between  the  keel 
and  the  wave  surface  (equation  (15) ) . 
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VII.   CONCLUSION 

The  work  presented  here  develops  a  slam  force  prediction 
method  which  is  sensitive  to  seakeeping  parameters  based  on 
a  model  of  slamming  different  from  other  methods  developed 
to  date.   The  results  obtained  show  slam  impulses  and  time 
durations  less  than  those  in  the  Ochi-Motter  method  [14] ,  but 
impulses  greater  than  the  100  ton-sec  value  used  by  Kline  [3] 
It  is  felt  that  values  closer  to  those  predicted  here  should 
be  used  in  the  overall  A.B.S.  study  on  hull  stiffness.   As  a 
minimum,  the  results  of  this  work  have  shown  that  the  100  ton- 
sec  value  used  in  the  work  to  date  is  in  question  and  should 
not  be  applied  to  all  ships  under  every  operating  condition. 
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APPENDIX  A  -  LOCAL  SLAM  PRESSURE,  DURATION  OF  PULSE 

Ochi  predicts  that  the  slam  pulse  is  triangular  and  the 
duration  is  0.1  seconds  for  a  520  foot  ship.   From  this  he 
Froude  sacles  and  concludes  that  equation  (13)  holds  -  i.e., 


t1   =    4.4  x  10  3  /L 


From  full  scale  measurements  on  the  USCGC  UNIMAC  [19] , 
L  =  300  feet,  values  of  slam  duration  locally  were  of  the 
order  of  0.020  to  0.025  seconds.   Using  equation  (13),  this 
length  ship  results  in  a  predicted  duration  of  t.  =  0.076 
seconds.   This  discrepancy  was  too  great  to  ignore  without 
further  consideration. 

From  drop  tests  done  by  Chuang  [8],  [9],  and  [20],  on 
flat  bottom  impact  with  zero  deadrise  angle,  the  following 
equation  holds  for  deadrise  angles  less  than  3° 


t,  =  4  1  /C  .  (22a) 

1      o   air 


where 


1   =  half-width  of  infinitely  long  flat  plate 
o 

C  .   =  speed  of  sound  in  air  =  1125  ft/sec 
air 


For  the  UNIMAC  and  the  durations  of  slams  measured  this 
would  predict  the  half-width  at  frame  23  to  be  of  the  order 
of  6  feet.   Lines  were  not  available  to  determine  how  close 
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this  agrees  to  the  half-breadth  of  the  1/10  draft  at  this 
section,  but  it  is  a  reasonable  value  for  a  ship  with  beam 
41  ft. 

In  other  drop  test  data  [21]  equation  (22a)  held  quite 
well  for  many  different  sized  test  specimens  if  the  beam  was 
used  as  the  scaling  parameter.   Without  more  information  from 
full  scale  tests,  equation  (22a)  cannot  be  verified  and  there- 
fore will  not  be  used.   It  is  felt  that  the  local  half-breadth 
should  be  the  scaling  parameter,  not  the  length  as  used  by 
Ochi. 

It  is  proposed  here  that 


t   =  3.6  x  10  3  B'  (14) 


■1 


which  comes  from  equation  (22a)  would  be  a  better  estimate 
of  slam  duration  with  possible  modification  of  the  constant 
from  full  scale  data. 

Equation  (22a)  holds  for  small  deadrise  angles  less  than 
3°  where  the  compressibility  of  air  influences  the  duration  of 
pressure  pulse.   For  larger  angles  the  impact  pressure  time 
history  becomes  more  steep  fronted  with  about  half  the  duration 
predicted  by  equations  (22a)  and  (14) .   Therefore  above  3° 
equation  (22b)  holds  fairly  well. 


t,  =  2  1  /C  .  (22b) 

1      o   air 
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Although  equations  (22a)  and  (22b)  are  for  deadrise  angles,  a 
similar  result  might  be  expected  for  angles  between  the  keel 
and  the  wave  surface.   Therefore  equation  (22a)  is  probably 
valid  in  the  flat  region  near  the  slam  center  while  equation 
(22b)  is  probably  valid  outside  this  region. 

For  the  MARINER  equation  (14)  yields  values  of  t,  =  0.022 
at  half-station  4  to  t,  =  0.115  at  half-station  16.   Therefore 
it  results  in  much  lower  impulses  than  would  result  using 
equation  (13) .   For  the  FOTINI-L  larger  values  were  found 
using  equation  (14)  than  from  equation  (13).   At  half-stations 
4  and  10,  t..  =  0.130  and  0.189  respectively  using  (14).   This 
would  result  in  larger  impulses  than  one  would  get  if  (13) 
were  used.   If  equation  (22b)  were  applied  in  the  appropriate 
regions,  the  above  resulting  impulses  could  be  reduced. 

Use  of  equations  (22a),  (22b),  and  (14)  are  not  justified 
at  this  time,  but  it  is  felt  that  full  scale  data  should  be 
analyzed  to  determine  if  equations  of  this  form  are  valid.   If 
they  are,  further  reduction  of  calculated  impulses  by  the 
program  is  possible  in  many  ship  cases. 
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APPENDIX  B  -  COMPUTER  PROGRAM  LISTING 

This  appendix  contains  a  listing  of  the  FORTRAN  computer 
program.   It  was  written  for  an  Interdata  Model  70  computer 
and  requires  approximately  60  K  of  memory.   Two  subroutines 
are  called  which  are  not  listed  here  but  are  readily  available 
to  most  users  as  part  of  the  IBM  Scientific  Subroutine  Package 
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C      PROGRAM  SLAMCW  WALKER      JANUARY  1976 

C      PROGRAM   MAIN  THIS  PROGRAM  CALCULATES  THE  SLAM  FORCES  ANO 

C      IMPULSES  GENERATEO  ON  A  SHIP  IN  AN  IRREGULAR  SEA 

DIMENSION  XIN(5),YIN(5) .XOUT(ll) .YOUT(ll) ,XX  ( 17. 1 1 1 . YY ( 1 1 ) .T (20) . 
C  VK(17).Y0(11).RM<17)  ,PV(17)  .EPPES(17).PRr>RPR(l7>, 

C  PRESU7.11)  .FOPCE(17.?0),C>VELtl7)  ,  TIMF  (17) 

DIMENSION  A (17) .OUR (17) .SLAMS (17) ,TM(20) .SHAPE (17) . 
CPRSLAMU7)  ,EMAX(17)  .  AIMPLS  ( 1  7)  ,PVEL  ( 1  7)  ,  XG  ( 1  7. 1 1 )  .XSCL  (4) 
OATA  ESL.ESP.NIN.N0UT.PI/99..R9..5.11.3.141592654/ 
DATA  N.M/8.S/ 

DATA  SHAPE /0. 89375, 0.791 6667, 0.69375. 0.6.0. 51 04167,0.425.0. 3475. 
CO. 266666, 0. 1 9375, 0.1 25, 0.1 6041 666. 0.00,-0. 05625, -0.108333. -0.15625 
C, -0.2,-0. 239583/ 
1   READ  (N.1000)  101, 102, 104, 108 

1000  FORMAT  (715) 

IF  (I02.LT.1)  GO  TO  999 
C      REAO  SHIP  NAME  AND  LOADING  CONOITION 
READ  (N,1001)  (T(I) ,1=1,20) 

1001  FORMAT  (20A4) 

C      READ  SHIP  DATA  CONSTANT  IN  PROBLEM 
READ  (N,1002)  RB,HH,ALBP«ALWL,STASP 

1002  FORMAT  (5F10.5) 
WRITE  (M.1003) 

1003  FORMAT  (M»,35X,»THIS  PROGRAM  CALCULATES  THE  SLAM  FORCE  DISTRIBUTI 
CON  VS  TIMEt,/,38X,«0N  THE  LOWER  1/10  OF  THE  SHIP  FROM  STATION  8  FO 
CRWARD.t ,////) 

C      WRITE  SHI°  NAME  AND  DATA 

WRITE  (MtlO04)  (T(I) ,1=1,20) ,BB,HH,ALBP,ALWL,STASP 

1004  FORMAT  (20X.20A4,//,30X.'H4LF  BEAM  =  • ,F  10.5, 10X . 'DESIGN  ORAFT  =  • 
C,*"10.5,//,30X.«LBP  =  «.F10. 2. 16X, "LENGTH  LOAO  WATER  LINE  =  '.F10.2 
C.//.30X, 'DISTANCE  BETWEEN  STATIONS  =  ' ,F6.2,40X,//,35X, • ( ALL  OtSTA 
CNCES  IN  FEET) •,////) 

C      CALCULATE  THE  K  VALUES  USING  SUBROUTINE  SPOO 
H=HH/10. 
S=HH/100. 
XOUT<1)=0.0 

00  10  1=2,11 
XOUT(I)=XOUT(I-l)  ♦  S 

10  YY(I)=XOUT(I) 

1  =  0 

C      READ  OFFSETS  AT  ANY  5  WATERLINES 

REAO  (N,1005)  <XIN( J) ,J=1,5) 
1005   FORMAT  (5F10.5) 

20   READ  (N.100S)  (YIN(J) ,J=1.5) 

1*1*1 

IF  (YIN(5).LE.0.5)  GO  TO  30 
C 

C      CALCULATE  OFFSETS  IN  10  INCREMENTS  FROM  B.L.  TO  1/10  DESIGN  DRAFT 
C      USING  SUBROUTINE  UGLYOK 


C 


CALL  UGLYDK  (NIN.NOUT.XIN, YIN»XOUT,YOUT.ESL»ESR) 
DO  21  J=l»ll 
JJ=12-J 
YO(J)=YOUT(JJ) 
21  XX(I,J)=YOUT(J) 
B=YOUT(ll) 
IAOP=I01 

CALCULATE  THE  FORM  COEF.  OF  THIS  STATION 
CALL  SPOO  (YO,VKK,B.H,S.BB.HH.IAOP.I) 
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30 

4003 
C 


1011 


VK(I)=VKK 

IF  (I04.LT.1)  GO  TO  20 
D=FL0AT(I)/2. 

WRITE  (M,1006>  0, <YIN< J) »XIN(J) 
C         XOUT(J) ,J=6.11 ) 

1006  FORMAT  (M»»<>9X.50(»»»),/»50X,» 
CF5.lt*  OFFSETS  (FEETM.  9X.»<>« 
C50X,»««,A«X.«»»»/»50X»«»,.9X.»  t 
C*/»S0Xt •••,?3X,»<»»,?<»X.  •••,/,  50 
C»llX.»<»».12X.«»«,/,5nx.«,,,5X.» 
C»«t5Xt»Y»  «6X. •••,/» SOX, 50 (•••) , 
C2X. •••./.S< SOX. •••.F8.2,3X,t»t, 
C./> .6(50X.i»»,llX.»»«.llX.»»».F 
C»> ./.'O') 

60  TO  20 
NLINFS=I-1 

REAO  tN.4003)  IPLT,  (XSCL(K),K= 

FORMAT  (I5.4F10.2) 

CALCULATE  AN  AVERAGE  FORM  COEF. 

CALL  KEFFEC  (VK.AVGK) 

WRITE  (M.1011)  AVGK 

FORMAT  (  »0»,39X»50(  •••)  ./.<»0X.' 
C  6X. 'AVERAGE  HULL  SHAPF  FACTOR 
C40X»,»« ,4<*x,  t», ./.40X.S0I •••)./ 

READ  IN  INTEGERS  FOR  CONTROL  OF 
100  REAO  (N.1000)  109, 103. 105, 106. I 

IF  (I09.LT. 1)  GO  TO  1 

REAO  (N.1007)  (T(I) .1=1,20) 

1007  FORMAT  (?na4) 
NSTA=NLINES 

WRITE  (M,1020)  (T(I) ,1=1,20) 
1020  FORMAT  (»1».1X.20A4.//) 

REAO  (N.1008)  FROUOE. SPEED. SEA. 

1008  FORMAT  (8F10.5) 
READ  (N.400P)  TSTORM 
FORMAT  (F10.5) 
HSIG=SEA»ALBP 

IF  (TSTORM. GT. 0.01)  GO  TO  501 
IF  (HSIG.LE.35.)  GO  TO  500 
TSTORM=S0.0-1.0«HSIG 
GO  TO  501 

TST0RM=78.n-l ,8«HSIG 
IF  (HSIG.GE.20.)  GO  TO  501 
TST0RM=45.-0.15»HSIG 

501  CONTINUF 

IF  (FROUOE. GT. 0.00001)  GO  TO  11 

FROUDE  =  SPFE0»1.6889/SORT(32.1 
110  SPEED  =  FR0UDE»S0RT(32.l74»ALWL 

FSCALE=SORT(ALPP/520.) 

THETA  =(HA-HF)/ALBP 

VAFTS=VAFT«FSCALE 

VTHRS=VTHR»FSCALE 

ALB=ALBP/(?.«B9> 

ABT=4.«BB/(HA*HF) 

IF  (IOS.EO.l)  GO  TO  103 

REAO  (N.1008)  (RM(I) .PV(I) .1=2. 

CALCULATE  THE  RELATIVE  MOTIONS 

CALL  RELMV  (RM.RV) 

GO  TO  10* 


.YOUT(J) .XOUT(J) ,J=1,5) ,(YOUT(J) , 

♦•,48X,»*»,/,50X»«»»,1 IX, 'STATION', 
«/, SO X.  ••».<»PX. ••',/, SOX, 50 (•••),/, 
NPUT' ,9X,'»',7X,'CaLCULATEO',7X, ••• 
X, 50 (•••)./. 50 X. •••, 1 1*. •••» 1 IX. ••• 
X«»5Xt»*««5X««Y».5X,»#».5X.,X»,5Xf« 
/,50X,,«f,llX,»«»,UX,«««,UX,«»«,l 
F8.?,3X,»»«,F8.2.3X,«»»,F8.2,4X,»»» 
8.2»3X»»»«,FB.2,4X,«»»»/)»50X.50(»» 


1,4) 
FOR  FORWARD  PORTION  OF  SHIP 


•i,48X,»»«,/,40X,»»«, 
=  • .F8.5.6X, ••»,/, 
/) 

OUTPET  OPTIONS 
07,1010 


HF,HA,CB,VAFT,VTHR 


4000 


500 


74»ALWL) 
)/1.6889 


8,6) 

AND  VELOCITIES  AT  ALL  HALF-STATIONS 
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103    NR=0 

122    IF(NR.GE.NSTA)  GO  TO  1 P4 

READ  (N.1008)  RM(NR*l).RV(NR*l),RM(NR*2),RV<NR*2),RM<NRO>, 
CRV(NRO)  ,RM(NR»4)  ,RV(NR*4) 

NRbNR«4 

GO  TO  122 
10*    CONTINUE 

C      CALCULATE  SLAM  PRESSURES  AT  THE  KEEL.  NUMBER  OF  Si AMS.  EXTENT  OF 
C      SLAM.  VELOCITIES  OF  SLAM.  AND  PROBABILITY  OF  SLAM 

CALL  SLAMPP  (RM.RV.HF.HA.ALWL.NSTA.EPRES.VK.SLAMS.PROBPR.OVEL.VTHR 
CS.TSTORM.HSIG.PRSLAM,PVFL.STASP> 

NCENTR  =  NSTA/2  »1 

IF     (VAFT.GT.0.01)    GO    TO    105 

IF    <NSTA.FO.0)    GO    TO    105 

Hl*(l.-SIN( (PI/2.)-0.071<»«FLOAT(NSTA>))»MSIG/2. 

IF     (IO6.EQ.0)    GO    TO    10<S 

TT«H1/PVEL (NCENTR) 

GO    TO    107 

106  TT=Hl/DVEL (NCENTR) 

107  VAFTS=FlOaT(NSTA)»STASP/(4.»TT) 
VAFT=VAFTS/FSCALE 

105    IF  (NSTA.EQ.0.ANO.IO3.EO.1)  READ  (N.1089)  JJJ 
1089   FORMAT  (13) 

WRITE  (M,im<0  CB.FROUOE. SPEED. AL8.ABT. SEA, HSIG.TSTORM.THFTA.HF, HA 
C.NLINES.NSTA 
1014  FORMAT  ('0',1<?X.'THE  SHIP  OPERATING  CONDITIONS  ARF:',//, 
C20X, 'BLOCK  COEF.  =  '.F7.4./. 

C20X,'FROUDE  NO.  =  « .F7.4.25X. 'SHIP  VELOCITY  =  ',F<S.2,»  KNOTS'./. 
C20X. 'LENGTH  TO  BEAM  PATIO  =  ',F6.3,/, 
C20X,'BEAM  TO  DRAFT  RATIO  =  '.FiS.3,/, 

C20X, 'NON-DIMENSIONAL  SEA  STATF  =  • .F7. 4t 10X. 'SIGNIFICANT  WAVE  HEIG 
CHT  s  '.F5.C  FEET'./. 

C20X,'TIME  OF  STORM  =  ',F5.0.'  HOURS'./. 
C20X.'TRIM  ANGLE  (RADIANS)  =  '.F7.4./. 

C20X, 'DRAFT  FORWARD  =  • .F5.2.24X . 'DRAFT  AFT  =  '.F5.2./, 
C20X.'NUMBER  OF  SHIPS  LINES  READ  =  '. 13. 13X. 'NUMBER  OF  HALF  STATION 
CS  CONSIOEPEO  =  '.13./) 
WRITE  (M.1015)  VAFT.VAFTS.VTHR.VTHRS 
1015   FORMAT  (20X. 'TRAVELING  VELOCITY  =  ',F9.1,»  FT/SEC  FOR  520  FT  SHIP 
C'»//. 

C20X. 'TRAVELING  VELOCITY  =  '.F9.1,'  FT/SEC  FOR  THIS  SHIP'//, 
C20X,'THPESHOL0  VELOCITY  =  '.F6.2.'  FT/SEC  FOR  520  FT  SHIP'.//, 
C20X, 'THRESHOLD  VELOCITY  =  '.F6.1.'  FT/SEC  FOR  THIS  SHIP',//, 
C20X,'THE  PROBABILITY  OF  EXCEEDING  THESE  SLAM  FORCFS  (ALPHA)  =  0.01 
C'.//) 
IF  (NSTA.EO.O)  WRITE  (M.3030) 
IF  (NSTA.FO.O)  GO  TO  100 
3030  FORMAT  COX.'NO  SLAM  CALCULATIONS  PERFORMED  SINCE  THE',/.40X, 
C'NUMBER  OF  SLAMS  IS  LESS  THAN  1  FOR  ALL  STATIONS',//) 
IF  (IO6.EQ.0)  GO  TO  150 
C      CALCULATE  SPATIAL  DISTRIBUTION  OF  PRESSURE  ALONG  GIRTH 
CALL  PRESS  (PROBPR,NSTA,PRES,XX,YY.S.XG) 
GO  TO  63 
150  CALL  PRESS  (EPRES.NSTA ,PRES,XX,YY,S,XG) 
63  CONTINUE 
C      INTEGRATE  PRESSURE  TO  GET  FORCES  FOR  EACH  HALF-STATION 

CALL  FORCES  ( FORCE. NST A. TIME .DELTAT.PRES.XG.ALBP.nvEL.VAFTS.STASP. 
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C  YY.PVEL.I06) 

C      CALCULATE  HALF-SINE  APPROXIMATION  WITH  EQUAL  IMPULSE  ANO  DURATION 
C      AT  EACH  HALT-STATION 

CALL  APPRO*  (  NSTA, FORCE .FMAX. AIMPLS. A.DUR.DELTAT) 

IF  (1010. EO. 2)  GO  TO  30?0 

IE  (1010. EO.l)  WRITE(M,3012) 
1012  FORMAT  (M». 126  <•••)./. IX. 2<"m,9X).'#'.14X.'»'.17X.'»'.6(11X.'»'> 
Ct/»IX.2('»',9X) .'•','   BEL  MOTION   •   EXT   PRESSURE   •♦.llXt 
C'»   MAXIMUM   ••tllXt'*  AMPLITUDE  •',11X,'<»     TIME    •'./. 
C1X."    HALF   •    FORM   «,34< •«•>,'    NUMBER   •     SLAM    •     SLA 
CM    •   OF  SINE   •     SLAM    •     SLAM    •»./. 

C1X.'*  STATION  •    COEF   •  REL  VELOCITY  •  MOST  PROB  PRESS  •   OF  SLA 
CMS  •    FORCF    •   IMPULSE   •    APPROX   •   DURATION  •    STARTS   ••. 
C/tlX.'»'.9X.'»    (K)    •    (NON-OIM)   •»,6X,« (PSI) • »6X, •• • . 1 IX  . 
C"    (KIPS)   •  (KIP-SEC)  •    (KIPS)   ••t2M    (SEC)    ••>./♦ 
C1X,126<"")) 

IF  (IO10.EO.1)  WRITE  (M.3013)  ( J.VK ( J) ,RM ( J) .EPRES ( J) .SLAMS ( J)  . 
CFMAX(J) .AIMPLS(J) .A (J) .OUR (J) .TIME < J) .RV ( J) .PROBPP(J) .J=1,NSTA) 

3013  F0RMAT(17(lx.'»',3X.I3,3X.'»',2X.F6.3.'  •• .4X.F6.3.4X, •  •' .6X.F6. 1 . 

csx. 

C'#'.2(1X.F9.0.1X,'»'> ,2X.F8.1,1X,'«M,2X.F8.0.1X.»»'.2(3X.F6.3.2X, 
C'#»> ./.IX,'«'.2(9X,'»») ,4X,F6.3.4X.'«''.6X.F6.1.5X.'»'.6(llX,'»'>. 
C/»1X»126('«'> ./) ) 

IF  (1010. EO.O)  GO  TO  3015 

IF  (106. EO.O)  WRITE  (M.3011) 

IF  (106. EO.l)  WRITE  (M.30U) 
3011  FORMAT  <30X.'»««  EXTRFME  PRESSURE  USED  IN  FORCE  CALCULATIONS  •••». 
C//) 

3014  FORMAT  <30x.'»«»  MOST  PROBABLE  PRESSURE  USED  IN  FORCE  CALCULATIONS 
C  ••••»//) 

GO  TO  3020 

3015  DO  3000  KM=1,NSTA 
TM(1)=TIMR(KM) 

DO  3006  KN=2,20 
3006  TM(KN)=  TM(KN-l)  ♦  DELTAT 

WRITE  (M.3001)  KM 

IF  (107. EO.O)  GO  TO  3010 

IF  (106. EO.O)  WRITE  (M.3002) 

IF  (106. EO.l)  WRITE  (M.3004) 

WRITE  (M.3005)  (TM(J) .FORCE(KM.J) .J=1.20) 

WRITE  (M.30O7) 
3010  CONTINUE 

WRITE  (MOOnfl)  EPRES(KM)  .PROBPR(KM)  ,VK(KM) 

WRITE  (M.3009)  RM(KM) .RV<KM> .SLAMS(KM) .PRSLAM(KM) .DVEL(KM) ♦ 
C       PVEL(KM) .FMAX(KM) .AlMPLS(KM) ,A(KM) .OUR (KM) .TIME (KM) 
3009  FORMAT  (30X, 'NON-DIMENSIONAL  RELATIVE  MOTION  =  • ,F7.4,/, 
C30X, 'NON-DIMENSIONAL  RELATIVE  VELOCITY  =  «,F7. 4,/, 
C30X, 'NUMBER  OF  SL*MS  DUPING  EXTREME  STORM  =  '.FP.n./. 
C30X, 'PROBABILITY  OF  SLAM  IMPACT  =  '.F8.5./. 
C30X,'EXTREMF  RELATIVE  VFLOCITY  =  ',E8.1.'  FT/SEC',/. 
C30X.'MOST  PROBABLE  RELATIVE  VELECITY  ■  '.F8.1.'  FT/SEC'.//. 
C30X, 'MAXIMUM  VALUE  OF  SLAM  FORCE  =  'tF10.lt*  KIPS'. /t 
C30X, 'IMPULSE  VALUE  =  '.FlO.O.'  KIP-SEC'./. 

C30X. 'AMPLITUDE  OF  SINE  APPROXIMATION  =  '.F10.1.'  KIPS'./. 
C30X, 'DURATION  OF  SLAM  IMPACT  FORCE  =  '.F8.5.'  SEC'./. 
C30X,'TIME  SLAM  STARTS  =  '.F8.5t«  SEC'.//) 

IF  (106. EO.O)  WRITE  (M.30U) 

IF  (106. EO.l)  WRITE  (M.3014) 
3000  CONTINUE 
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3020  CONTINUE 

3001  FORMAT  ('l'.29x.40('»') ,/.30X.»»», 38X. »••»/, 30X. ••• . 10X. 'HALF  STAT 
CION   •.I3.11X.""./,30x.'°'.3flX.'»'  , /.30X.40('»« ) > 

3002  FORMAT  (  30X  ,  •  •• .  1  7X.  ••  •  .3X  .  'EXTREME  FORCE  •  ,3X  .  •  ••  ./.30X.  •••  .3X.  «T 
CIHE  (SEC) ».3X,'<»',5X,'F0R  DESIGN' .SX . ••' ,/,30X. ••'. 17Xt ••• .7X, • (K 
CIPS) '.7X.'o'./.30X, 40 <•••)./. 30 X,'«".17X.'»'. 20 X.*«») 

300*  FORMAT  (30x«'»'.17X»'»'.3X,'MOST  PROBABLE' ,3X» '•. ,/,3ox, ••«. 3X, • T 
CIME  (SEC) • »3X.'»'.5X,«  FORCE  • .SXi ••• »/,30X • •  •» t 17Xt ••• ,7X. • (K 
CIPS)'.7X. •»',/, 30 X, 40 (•*•>, /«30X.'»»»17X.'»». 20 X»«»»> 

3005  FORMAT  (30X»»««.5X,F7.3.5X,«<".*X,F8.1»6X«««»t) 

3007  FORMAT  <30X,&0( '»•> .//) 

3008  FORMAT  (30X. 'EXTREME  PRESSURE  =  '.F10.0.'  PSI •./. SOX. 'MOST  PROBABL 
CE  PRESSURE  =  '.F10.0.'  PSI • »/«30X, 'FORM  COEF.  <K>  s'.FlO.S./) 

TIMPLS=0.0 
DO  200  I=1.NSTA 
200   TIMPLS=TIMPLS  ♦  AIMPLS(T> 

TSLAMT=TlME(NCENTR)*OUR<NCENTR> 

TTTsTIMPLS 

TA  «=  PI»TIMPLS/(2.«TSLAMT) 

L0C1=NCENTR 

WRITE  (M.1017)  TA.TSLAMT 

1017  FORMAT  (20X.»  AVERAGE  SLAM  FORCE  =  • »F8.0. '"SINE <PI#T/» t F6.4. • ) 
C  KIPS',//) 

WRITE  (M.101«)  LOCI 

1018  FORMAT  (25X,'MI0P0INT  LOCATION  IS  HALF-STATION' .14.///) 
TIMPLS  =  0.0 

00  202  I=1.NSTA 
202  TIMPLS=TIMPLS*AIMPLS(I)«SHAPE(I) 
TB=RI»TIMPLS/(2.«TSLAMT) 
WRITE  (M. 10191  TB.TSLAMT 

1019  FORMAT  (20X, 'MODAL  AVERAGE  SLAM  FORCE  =  • »F8.0 • ''SINE (PI»T/« .F6.4 
C»»)      KIPS'//) 

WRITE  (M.1018)  LOCI 

WRITE  (M.1069)  TTT 
1069  FORMAT  ( '0 ', 10X . 'TOTAL  SLAM  IMPULSE  =  '.F8.0.'  KIP-SEC »///) 

IF  (IO3.EQ.0)  GO  TO  100 
C      CALL  PLOTTING  ROUTINE  FOR  TIME  HISTORY  OF  FORCES 

CALL  PLOT  (TIME. FORCE. OELTAT.NSTA.IPLT.XSCL) 

GO  TO  100 
999  WRITE  (M.1010) 
1010  FORMAT  CO'.'  END  OF  JOPM 

STOP 

END 
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SUBROUTINE  UGLYDK  (NIN.NOUT .XIN. YIN* XOUT t YOUT tESL .ESR) 
C      THIS  SUBROUTINE  TAKES  THE  OFFSETS  AT  EACH  INPUT  WATERLTNE  AND 
C      USES  SPLINE-CUBIC  INTERPOLATION  TO  DETEPMlNE  THE  10  EOUALl Y 
C      SPACED  WATE«LINE  OFFSETS  BETWEEN  RL  AND  1/10  DESIGN  DPAFT 

REAL#B  H.D«A,S.AF»BE«CE«DE»HALFtTwO«SIXtPAO,SLP.DSIN,DCOS 

DIMENSION  <IN(S).YIN(S)  ,XOUT(l 1)  .YOUT(ll) tH(5)tO<S).A(25>.S<5). 
C  AE (5) ,BE(5) .CE<5) »DE(5) 

DATA  HALF/0. 5D00/.T WO/2. 0D00/.S I X/6.0D00/. RAD/ 1.745329250-02/ 

NH1=NIN-1 

NM2=NM1-1 

NE0=NH2 

DO  1  N=l.NMl 

H(N)=XIN(N*1 )-XlN(N) 

1  0(N)=(YIN(N»1 )-YlN(N) ) /H(N) 
IF  IESL.GT.90.0)  GO  TO  2 
NEO=NEQ*l 

2  IF  (ESR. GE. 90.0)  GO  TO  3 
NEQ=NEO»l 

3  NSO  =  NEO«<>2 
DO  4  N=l .NSO 

4  A<N)=0.0 
J=l 

1  =  1 

IF  (ESL.GE.90.0)  GO  TO  ft 

A(l)=TWO°H(l ) 

A(2)=H(1) 

SLP=ESL»RAD 

S<1  )  =  (D(1)-(DSIN(SLP)/DC0S(SLP)  >  >»SIX 

J  =  J»1 

M=L*NEO 

A(M)=H(1 ) 

L=M»1 

6  DO  5  N=l ,N*2 

IF  (N.GT.l)  A(L-D=H(N) 

A(L)=TWO#(H(N)  ♦HCNM)  ) 

IF  (N.LT.NH?)  A(L«1)=H(N*1) 

S<J)  =  <D(NM)-D<N)  )»SIX 

J=J«1 

5  L=L»NEO*l 

IF  (ESP. GE. 90.0)  GO  TO  7 

A(L-1)=H(NM1) 

A  <L)=-TW0»M(NM1 ) 

L=L-NEO 

A(L)=-H(NM1 ) 

SLP=ESR»RAD 

S(J)=(D(NM1)» (OSIN(SLP) /OCOS(SLP) ) )*SIX 

7  CALL  SIMQ  (A,S.NEO»KEPROR) 
IF  (ESL.LT.90.)  GO  TO  °t 
DO  9  N=] .N"2 

M=NM2-N*2 

9  S(M)=S(M-1) 

S(1)=0.0 

8  IF  (ESR.GT.90.0)  S(NIN)=0.0 
DO  10  N=ItNMl 

AE(N)=(S<N*1)-S(N))/(SIX»H(M)) 
BE (N)=HALF«S(N) 

CE (N)=D(N)-M(N)«( TWO'S (N)*S(N*1) )/SIX 
10  DE(N)=YIN(N) 
J=l 
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DO  11  N=l.NOUT 

IF  (XOUT(N) .GE.XIN<2> )  GO  TO  14 

J=l 

GO  TO  15 
\i*    IF     (XOUT(N)  .LT.XIN(NMl)  )  GO  TO  16 

J=NMl 

GO  TO  15 
16   IF  (XOUT(N) ,GF.XIN(J*1 > )  GO  TO  12 
15  Hl=XOUT (N)-XIN(J) 

H2=Hl»»2 

H3=H1»H2 

YOUT(N)=AE< J) •h3»BF(J)«h2*CE(J)»H1*DE(J) 

GO  TO  11 
12  J=J*1 

GO  TO  16 
11  CONTINUE 

RETURN 

END 
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SUBROUTINE  SPOO  (X . VKK.B.H, S.BB.HH. I AOP, ISTA) 

SUBROUTINE  SPOO  ( X , VKK .R) H,S,BB,HH, I AOP, ISTA) 
C      THIS  PROGRAM  CALCULATES  THE  FOR"  COE*".  K  DEVELOPED  BY  OCHT  TO  BE 
C      USED  IN  PREDICTING  THF  SLAM  CHARACTERISTICS  OF"  A  SHIP  FROM  THE 
C      SHIPS  LINES.      REFERENCE  SNAME  TRANSACTIONS  1R73 

DIMENSION         X(]l )  tXM(ll) ♦YM(11 ) , ZM  C 1 1 ) .ROOT (5)  .W(f>)  . 
2  XN(ll)  ,  7N  ( 1 1 )  .XI  (11)  ,X2(11  )  ,Z1  (11)  .Z2M1> 

2  .WW (6) ,ROOTR(S> .ROOT! (5)  .C0L(M 

DATA  P I /3.1M5R2f>/,  EPS/0.  0005/,  N/R/.H/5/ 

IF  (IAOP.EO.1)  WPITE(M.lOOl)  B.H,B9,HH 

IF  (IAOP.FO.1)  WRITE(M.lOll)  (X<I)  .1  =  1,11) 

IF  (IAOP.EO.1)  WPITF  (M.1002)  ISTA 
1002   FORMAT  (28X,»  HALF-STATION  NUMBER   ».I3./) 
C 
C      CALCULATE  AREA  AND  MOMENT  OF  INERTIA  BY  SIMPSON'S  RULE. 

C 

A  =  0. 

DO  10  1=2, «. 2 
10  A=A«^.«X(I) «2.»X(I*1) 

A=2.#S/3.#(A»X<1)  »X(U)«<».*XU0)> 

TI=0.0 

AI=0.0 

DO  20  1=2. «, 2 

AI=AI-1. 

TI=TI*a.#X(I)»(AI*S)»«2  ♦  2.»X(I*1)»((AI»1.)»S)*#? 
20  A1=AI»1. 

TI=  2.»S/3.»(TI*4.»X(10)»(9.«S)«»2*X(11)#(10.«S)»»2) 

C 

C      CETEPM!NE  M.A3.A5.  *ND  U. 
C 

DL=H/B 

IF  (ABS(DL-1.)  .LE.  EPS)  DL=DL*EPS 

A  =  A<»2./PI 

C2=TI«B./PI 

TK1=B»«2»(1 .♦DL*DL»«2) 

TK2=B»( 1 .-DL) 

TK3  =  -3.»B<» (1 .«DL) 

TK4=(A»TK1) /TK2 

TK5=TK3/TK? 

TK6=2./TK2 

Dl=DL»»2»A.»DL-7. 

D2=-DL«»3*DL°»2  -3.»DL*3. 

D3=-DL*5. 

D4=DL*#3^DL*»2*3.,»DLO. 

D5=-DL»»2-R.»DL-12. 

D6=3.»DL»7. 

D7  =  -l  .75«DL*»^-0.S»DL,*»3-DL»»2-1  .S'DL'-l^S 

W(1)=-4.»T<6 

W(2) =2.«B«TK6«D3-^.«TK5-6. 
W(3)=B«(2.»<D6*T,'5»D3)^0*TK6»Dl)-^.oTKft 
y  (4)=B«(2.»TK'<,«D3*B»(D5*TK5»D1*0.5»B#TK6*02)  ) 
W(5)=B»«2*(T'<^»Dl*0.S»o«<3.*D<»*TK5*D2)) 
V<6)  =0.5»B«*3»(T'«'4»D2*H»37)  -C2 
IF  (IAOP.EO.1)  WRITE  (M.3000)   (W(I). 1=1.6) 
3000  FORMAT  (•   W=  • .6 (E 15.5. 3X ) ♦ /) 
MM  =  5 

DO  30  1=1.* 
J=7-I 
WW(D=W(J) 
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30   CONTINUE 

CALL  BAIRS  (WW.COL.MM.ROOTR.ROOTI.KK) 
ERROR=0. 000001 
NR  =  0 

DO  40  1=1.5 
ROOT(I)=0.0 

IF  <A8S(R00TI  (I)  )  .GT. ERROR)  GO  TO  40 
IF  (ABS(ROOTRd)  )  .LT.FRR3R)  GO  TO  40 
*jR=NR»l 

ROOT (NR) =ROOTR( I) 
40  CONTINUE 

IF  (NR.EO.l)  GO  TO  39 

00  38  J=l,3 

KB=J»1 

DO  37  I=KB.NR 

IF  (ROOT(J)  .GT.ROOTU)  )  GO  TO  37 

R=ROOT(J) 

ROOT (J) =ROOT ( I) 

ROOT(I)=R 

37  CONTINUE 

38  CONTINUE 

IF  (NR.EQ.M  GO  TO  39 

IF  <R00T(4) .GT.R00T(5) )  GO  TO  39 

R=ROOT(4) 

R00T(4)=R00T (5) 

ROOT (5)=R 

39  CONTINUE 

IF  (IAOP.EO.I)  WRITE (M, 3001)  (ROOT ( I ) . I = 1 .5) .NR.KK 

3001  FORMAT  (•   ROOT  =  •«5(ri0.5t3X) t/t*  NR  =   «.I5.»   KK  =  «.l5./> 
IF  (KK.GE. 1 )  GO  TO  99P 

00  25  I=1,\P 

U=R0OT(I) 

ALFA=2.«U/P 

A5=TK4/U*TK5»TK6»U 

Al=  (1 .-DL) /ALFA-A5 

A3=(1.»0L)/ALFA-1. 

IF  (IAOP.EO.I)  WRITE  (",3002)  Al.A3.A5 

3002  FOPMAT  (•   Al  =  •  ,F10.5.3X,«A3  =  • .FlO .5.3X. • A5  =  «.F10.5,/) 
IF  (A3  .LT.  0.)  GO  TO  ?5 

IF  (ABS(Al-FPS)  .LE.  1.  .AND.  ABSU3-EPS)  .LE.  1.  .AND. 
2   ABS(A5-EPS)  .LE.  1.)  GO  TO  26 

25  CONTINUE 
GO  TO  998 

C 

C      COMPUTE  K-VALUE 

C 

26  TK=-3. 59893  ♦  2.4l988»Al  -  0.872855»A3  ♦  9.62395»A5 
VK=EXP(TK) 

VKK=VK 

IF(IAOP.GT.O)  WRITE  (M,?000)  Al . A3. A5.U. VK 
IF  (IAOP.GT.O)  WRITE  (M.2004) 
A=A*PI/4. 

IF  (IAOP.EO.I)  WRITE  (M.3005)  A 
3005   FORMAT  (1X,«   HALF-SECTION  AREA  OF  1/10  DRAFT  =  «,F10.5./) 
IF  (VK.LE.0.19)  GO  TO  109 
VKK=  0.027»8»B/A 
IF  (IAOP.EO.I)  WRITE  (M.3004)  VKK 
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3004  FORMAT  (lx,»  THE  REGRESSION  ANALYSIS  <  NOT  USEO  »t/tlX.»  K  =  • 

C»E10.6OX,i  WHERE  K  =  0.027  »B»8/A»./> 
109   CONTINUE 
60  TO  100 
998  WRITE  (M.2002) 
100  CONTINUE 
RETURN 
1001   FORMAT  (1H1,5BX,1?H«««  SPOO  •••   //28X .72HTHIS  PROGRAM  COMPUTES  T 
2HE  CONFORMflL  MAPPING  COFFFICIENTS  AND  THE  K-VALUE/28X, 
2  63HNEE0ED  TO  COMPUTE  SI.AMMING  PQF.SSJRE  BY  THE  RECESSION  EOUATION 
2.  //2ftX,l^HHALF  3PFAOTH  =  .F10.5.5H  FCcT,3X, 

2  25HHEIGHT  ABOVE  THE  BOTTOM  =  ,F10.5»5H  FEET « /28X , 1 4HHALF  B^AM 
2=  .F10.5»5h  FEET.21X,7Mn9flFT  =  .F10.5.5H  F^FT) 

1011   FORMAT  (/  2RX.31HTHE  VALUES  X  .X  X    ARE  .F 1 0 .  S/<»0X .  1H0 . 

2  3X,1H] ,7x,?Hl0.5X.  9(F10.5/5RX> ,F10.5.5H  FFET) 
2000   FORMAT  <///?Rx,74HTHF  COEFFICIENTS  OF    THE  TRANSFORMATION  OF  THE  C 
2IRCLE  IN  THE  ZFTA-PLANE  AR£//?8X.SHA   =  . E 1 5. 8. ?X ,SHA   =  .E15.8, 
2  7X.5HA   =  .E1S.«/29X, lHl .26x«lH3f26X»lH5///2«X, lQHTHr  SCALE  FACTO 
2R  IS.4lX»l4HTHr  K-VALUE  IS//28X.4WU  =   .E  15.8.  3<SX  ,  CHK  =   .E15.A) 
2002   FORMAT  (///?8x,*fHA  SOLUTION  FOR  U  HAS  NOT  BEEN  FOUND.   t>LFASE  CH 

2ECK  THE  INPUT  DATA.//) 
2004   FORMAT  ( /3<SX , S7HWHERF   <  =  EXR  (-  3.599  ♦  2.420  A   -  0.873  A   ♦  9 
2.624  A  )   /69X.1H1.10X.1H3.10X.1H5) 
END 
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SUBROUTINE  SLAMPR  (RM.RV.HF  ,HA  ,  ALWL  ,  ISTA  ,EPPES,  VK  ,SL  AMS  ,PPORPR, 
C  DVEL»VTHRS,TSTq:?m,HSIG,P}SlAM,PvFL»STaSP) 

C      THIS  SUBROUTINE  CALCULATES  THE  EXTREME  PRESSURE  AT  EACH  STATION 
C      WHERE  THE  PROPARILITY  OF  SLAM  IMPACT  IS  GREATFR  T«AN  1 

01 MENS I  ON  PM(17)  ,PV(17) .EPRESU7) .VK<17> ♦PR0BPR(17) ,DVELU7), 
CPRSLAM(17)  .SLAMS (17)  ,PvrLU7) , X<17) 

DATA  ALPHA,  r-,p  1/0.01,32.  174,  3.  14 15926/ 

DATA  M,N/5,R/ 

RSTARS=VTHPS»»? 

AL=0.983»ALWL 

T=TSTORM 

NSTA=ISTA 

DO  10  I=1,MSTA 

H=( <HA-HF) /40.) »FLOAT(I) »HF 

DRM  =  2.«M <AL«RM(I)»0.fl3)»*2) 

0RV  =  2.»G»AL«MRV<  I >  »0.R3) »»2 

PRSLAMd  )  rrxPj-HoH/ORM-RSTARS/DRV) 

SLAMS ( I  J  =  <  3*00 . *T/ ( 2. »PI >  > »SQRT ( DRV/DRM) 'PPSLAM ( I ) 

IF  (SLAMS(I) .GE.l.)  GO  TO  11 

1STA=I-1 

NSTA=ISTA  ♦  1 

WRITE  (M.1002)  NSTA.TSTO^M.SLAMSU)  .HSIG 
1002  FORMAT  (•0«.»  THF  NUMRPR  OF  SLAMS  AT  HALF  STATION  «»I3,»  OURING  A 

cstopm  of  duration  »»F3.o««  hours  =  »«F7.4,/,i  fop  a  significant  w 

CAVE  HEIGHT  OF  »,F5.0.»  FEET',//) 

GO  TO  100 
11  EX=1./SLAMSU) 

EPRFS ( I ) =V<  (  I ) •  (RSTARS-n^V^ALOGC 1 .-( 1 .-ALPHA) ••FX) ) 

pRORPP(  T  )  =vx-  (I)»(RSTARS  ♦  ORV»ALOG(SLAMS(  I)  )  ) 

DVEL(I)=SO»T (EDRES(I)/VK(I) ) 

PVEL<I)=SQOT<PRO«PP<I)/V<<I) ) 
10  CONTINUE 

ISTA=NSTA 

WRITE  (M.lonl)  ISTA, SL AMS(ISTA) .TSTO^M, HSIG 
1001  FORMAT  (M'.'  ADDITIONAL  STATIONS  SHOULD  BE  CONSIOEPED  SlMCE  AT  HA 
CLF  STATION  «,I3,/,«  THF  DUMBER  OF  SLAMS  =  »,F5.1,t  FOP  A  STORM  OF 
CDURATION  *«F5.0f*  HOURS'*/.'  THE  SIGNIFICANT  WAVE  HEIGHT  =  »,F5.0t 
C  FEET',//) 

RETURN 
100    CONTINUE 

NCENTR=IST4/2*1 

HSTASP=STASP/2. 

XCENTR=FLOAT (NCENTR) «HST ASP-HST ASP/2. 

DO  1  J=1.NCFNTR 

X(J)=XCENTR- (FLOAT(J)»HSTASP-HSTASP/2.) 

1  CONTINUE 
J=NCENTR»1 

00  2  I=J,NSTA 
X(I)=X(I-1) «HSTASP 

2  CONTINUE 

DO  5  1=1 .ISTA 
A  =  -112.«HSIG<»X  (I)/(AL»<>2> 
EPRF.S(I)=EPPES(I)*(10.«»A) 
PROBPR( I )=DRORPR( I) »(10.#*A) 
5      CONTINUE 
RETURN 
END 
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SUBROUTINE  PELMV  (RM.RV) 
C      THIS  SUBROUTINE  TAKES  THE  INPUTEO  RELATIVE  MOTIONS  AND  VELOCITIES 
C      AT  STATIONS  1  K.  U    AND  ASSUMES  LINEAR  VARIATION  TO  DETERMINE  VALUES 
C       AT  THE  17  HALE-STATIONS 

DIMENSION  PM (17) .RV(17) 

SLOPE=(RM(B)-RM(?) )/6. 

B=RM(2)-SLOPE»?. 

DO  10  1  =  1 »  17 

RM(I)=SLOPE«fLOAT  <I)*B 
10  CONTINUE 

SLOPE  =  (PV(8)-RV(?) )/ft. 

B=RV(?)-SLOPE»?. 

DO  20    1=1.17 

RV(I)=SLOPF*FLOAT(I)  ♦  R 
20   CONTINUE 

RETURN 

END 
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SUBROUTINE  FORCES  (FORCE.NSTA.TIME.OELTAT.PRES.XX, ALBP.OVEL.VAFTS. 
C  STAS^.YY.PVEL.IOiS) 

C      THIS  SUBROUTINE  OOES  THF  SPATIAL  AMD  TEMPORAL  INTEGRATIONS  AT  EACH 
C      HALF-STATION  OF  THE  PRESSURES  TO  FIND  THE  SLAM  FORCE  TIME  HISTORY 

DIMENSION  FORCE (17,20) .PRES ( 17, 1 1 ) ,XX ( 17. 1 1 > .TIME ( 17) ,0VEL < 1 7) • 
C        YY(17)  .PVFLU7)  ,r(10)  ,FOR(17,20) 

TSLAM=0.1»SORT(ALBP/5?0.) 

OELTAT=TSLAM/10. 

N0ELT=INT(STASP/(2.»VAFTS»0ELTAT)) 

IF  (NOELT.LT.l)  NOELT=l 

T»STASP/(VAFTS«2.) 

NCENTR=NSTA/2*1 

TIME(1)=0.0 

00  A  I=2,NCENTR 

TIME(I)=TIME(I-1)*T 

4  CONTINUE 
I*NCENTR*1 

00  5  J=I.NSTA 
TIME(J)=TIME(J-1)-T 

5  CONTINUE 

00  20  I=1,NSTA 

IF  (I06.E0.1)  GO  TO  50 

OELY=OELTAT»OVEL(I) 

GO  TO  51 

50  OELY=OELTAT»PVEL(I) 

51  FORCE  (I,1)=0.0 
DO  21  J=2,5 

21  F0RCE(I,J)=(FL0AT(J-1)/5.>»(PRES(I,11»XX(I,1)  ♦  PRES<  I  .2)  MXXC I  .2) 
C-XX(I.l) ) ) 

00  22  J=6,10 

22  F0RCE(I,J)=(FL0AT(11-J)/5.)»(PRES(I,1)»XX(I,1)  ♦  PRES ( I .2) •UX ( 1.2 

O-XX  (I  .1  )  )) 

00  28  J=ll,20 
28  FORCE(I.J)=0.0 

00  23  J=3.11  t- 

F(1)=0.0  1 

OELX^tXXd.JJ-XXtl.J-lD^PRESd.J) 

00  2<»  JJ=2,5 
2*  F(JJ)=(FL0AT(JJ-1)/5.)»0ELX 

00  25  JJ=6,10  .:'   .  .    '« 

25  F  (JJ)  =  (FLOAT  (11-JJ)/5.)»0ELX  .  "    .  ;'.  .-• 
00  30  L=l,20 

M»INT ( YY ( J-l > / (DVEL ( I ) 'FLOAT (L) »OELT AT) ) 

INC*L 

IF  (M.LE.l)  GO  TO  31 

30  CONTINUE 
GO  TO  20 

31  M*INC 
MHsH.9 

IF    (MM.GT.20)    GO    TO    100 
DO   26    JJM.10 
K*M»JJ-1 

26  FORCE(I.K)=FORCE(I.K)     ♦    F(JJ) 
GO   TO   23 

100    IF    (M.GE.20)    GO    TO   20 
MM=20 
K«0 
00   27   JJ=M,MM 
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27  FORCE<I.JJ)=FORCE(I.JJ)*F(K) 
23  CONTINUE 

00  29  J=l,20 

FORCE<I.JJ=(STASP«FORCE(I.J)/1000.)»1<>4./FLOAT(NOELT> 
29    FOR(I.J)=FORCE(I«J) 
20  CONTINUE 

IF  (NDELT.LE.l)  RETURN 

00  60  I=1.NST4 

NEXT=1 

LAST*0 

00  70  M=2,N0ELT 

NEXT=NEXT*1 

LAST=LAST»1 

00  80  J=NExT,20 

JJ=J-LAST 

FORCE ( I . J) =FORCE (I.J) *FOR ( I . JJ) 
80   CONTINUE 
70   CONTINUE 
60   CONTINUE 

RETURN 

END 
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SUBROUTINE  PRESS  (EPRES.NSTA .PRES. XX ♦ YY.S. XG> 
C 

C      THIS  SUBROUTINE  DETERMINES  THE  SPATIAL  DISTRIBUTION  OF  THE  VERTICAL 
C      COMPONENT  OF  THE  PRESSURE  AS  A  FUNCTION  OF  DISTANCE  FROM  THE  CENTE»LINE 
C 

DIMENSION  EPRESU7)  .PRES<17tll)  .XX < 1 7. 1 1 >  * YY < 1 1 >  .CTHETA < l?)  ,X<11) 
C.XGC17.11) 

DO  10  i  =  1.n<;ta 

PRESU.l  )=EPRES(T) 

X(1)=XX(I,1 ) 

DO  12  J=1.10 

DELTAX=XX(T .J»1)-XX(I.J) 

HYP=SORT (S»»2*nELTAX»»?) 

X(J»1 )=X< J) »HYP 
1?   CTHETA( J) =OELTAX/HYP 

DO  13  J=ltll 
13  XG(I.J)=X(J) 

DELTAP=EPRES(I)/10. 

DO  1<»  J=2,ll 

A=FLOAT( J-l) 
1*  PRES(I»J)  =  (EPRES(I)-A«0ELTAO)»CTHETA(  J-l> 
10  CONTINUE 

RETURN 

END 
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SUBROUTINE  KEFFEC  (VK.AVGK) 
C      THIS  SUBROUTINE  CALCULATES  THE  AVERAGE  OF  THE  FIRST  11  HALF- 
C      STATIONS  CALLED  THE  AVERAGE  FORM  COEr. 

DIMENSION  VK<17) 

AVG*=0.0 

DO  20  J=l.ll 
20     AVGK=AVGK*VK(J) 

AVGK  =  AVGK/H. 

RETURN 

END 


SUBROUTINE  APPPOX  (NSTA .FORCE .FMAX . A  I M?LS» A .DUR.DFLTAT) 
C       THIS  SUBROUTINE  CALCULATES  FOR  EACH  HALF-STATION  THE  IMPULSE  BY 
C       INTEGRATING  THF  PORCE  TIME  HISTORY  AND  COMPUTES  A  HALF-SINE 
C       APPPOXIMATIONI  WITH  TH^  e;AME  IMPULSE  AND  DURATION 

DIMENSION  FO»CE(17.?0)  .FMAX  (17)  »ATMP;_S(17>  ,  M17)  fnuR(17>  .v(20) 

DATA  PI/3. 141592^.54/ 

DO  10  1=1. NSTA 

DO  11  J=10.?0 

IF  (FORCEM  .J)  .GE.O.OS)  GO  TO  11 

OUR ( I >  =FLOAT ( J- 1 ) *DELTAT 

NPTS=J 

GO  TO  12 
11  CONTINUE 

DUK(  I) =20. 'DEL  TAT 

12  00  13  J=1.?0 
Y(J)=F0RCE (T.J) 

13  CONTINUE 

CALL  MAXFOP  (Y.AMAX) 
FMAX (I ) =AMAX 
CALL  OSF  (OELTAT.Y.Y.20) 
AIMPLS(I) =Y(2P) 

A(I)  =  PI«»AIMPLS(I)/(2.»DUR(I)  ) 
10  CONTINUE 
RETURN 
END 
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SUBROUTINE  MAXFOP  (Y.AMAX) 
C      THIS  SUBROUTINE  DETERMINES  THE  MAXIMJM  FORCE  FROM  THE  FORCE  TIME 
C       HISTORY  AT  EACH  HALF-STATION  USING  A  PARABOLIC  FIT  THROUGH  THE 
C      HIGHEST  3  POINTS. 

DIMENSION  Y<20> 

FLAST  =  Y(?) 

DO  10  1=3. « 

N=I 

IF  (Y(I).LT. FLAST)  GO  TO  11 

FLAST  =  Y(I) 
10    CONTINUE 

AMAX  =  1.0E10 
11  AA=Y(N-2> 

BB=Y(N-1> 

CC=Y(N) 

B=0.5«(4.«PR-3.<»AA-CC) 

C=AA 

A=BB-AA-B 

X=-B/<?.«A) 

AMAX=A*X»*?  ♦  R»X  ♦  C 

RETURN 

END 


SUBROUTINE  PLOT  ( T IME . FORCES , DEL T AT , NST A . I »L T , XSCl  ) 
C      THIS  SUnWOuTINF  USES  THF  M.I.T.   J.C.P.    PICTR  PtOTTTNG  ROUTINE 
C      -  IF  THE  PROGRAM  Is  TO  fl[  USED  ON  ANOTHER  MACHINE.  THIS  IS  Tmf 
C      ONLY  SUBROUTINE  THAT  REO'JIRES  CHANGING  ALONG  WITH  THE  RELATEO 
C      CARDS  IN  main 

INTEGER*?  IXPLOT ( 17) .XL«3(40) 

DIMENSION  PL TR( 34.20) .TIME (17) .FORCES* 17.20) . xSCL(4> 

DATA  I XPLOT/1R, 19.20.21 . 22 .23.24 . 25. 2b. 27 .28,29.30 . 31 . 32. 33. 34/ 

DO  10  1=1. NST A 

T=TIME(I) 

DO  20  J=1.20 

PLTRU,  J)=FORCES(I,  J) 

11=1*17 

PLTRdI,  J)=T 
20   T=T»DELTAT 
10    CONTINUE 

CALL  PICTR  (PLTR. 34. XLAB. XSCL. NSTA. 20.-1.-1. 3. IPLT.0.1.IXPLOT) 

PAUSE 

RETURN 

END 
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SUBROUTINE    BA IRS ( ACOF , A , WRO,ROOTR,ROOTI , IER) 

REAL    ACOF() ) ,A(1) ,ROOTR(l) , ROOT  I (1) 

DATA    T0L2/1 .0E-3/.TOLX/S.0E-6/ 

DATA    TOLtlTMIN.ITMAX/l.OE-B.lO.lOO/ 

DATA    TOL3/1.0E-10/ 

N=NORO 

IF<N.LT.0)N=-N 

D=ACOF(NM  ) 

DOM=- 100.0 

IR  =  0 

C BASIC  ERROR  CHECKING 

IF(N) 1,1,2 
1      IR=1 

GO  TO  200 
?  IF(N-99)<»,4,3 

3  IR=2 

GO  TO  200 

4  IF(0)6.5,6 

5  IR=4 

GO  TO  200 

6  DO  10  1=1, N 

10    A(I)=ACOF(I)/0 
A(N*1 )=1 .0 

C REMOVE  TRIVIAL  ZEROES  FIRST 

J=2 

00  14  I  =  1.N 
IF(A(D)  12,11,12 
U    ROOTR(J-1)=0.0 
ROOTI ( J-i ) =0.0 
14    J=J*1 
C FIND  A  REAL  ROOT  IF  REDUCED  POLYNOM  IAL  IS  ODD 

12  LEFT=N-JM 
IF(LEFT) 100.90,13 

13  NN=N»1 
IJ=J-1 

C ESTIMATE  THE  MODULUS  OF  THE  LARGEST  ROOT 

DO  151  I=IJ.NN 
IF(AU)  .NF.0.0) 
100M=AMAXl (Al 0G(A9S(4<T)/A(IJ) ) ) /FLOAT ( \-\J* 1 ) ,DOM) 
151   CONTINUE 

DOM=EXP(DO*> 
OJ=l .0/DOM/DOM 
PJ=-2.0/DO* 
C       SCALE  THE  POLYNOMIAL  SO  LARGEST  ROOT  IS  NEAR  1.0 
DO  9  I=IJ,N 
DO  9  K=IJ,I 
9     A(K) =A<K) ZOOM 

IF(LEFT.NE.LEFT/2»2)GO  TO  15 
300   X0=0.0 

P0=A(J-1) 

C DETERMINE  IF  REAL  ROOT  IS  POSITIVE  OR  NEGATIVE 

Xl=1.0 

IF (P0.GE.0.0)X1=-1.0 

IP=1 

C EXPAND  RANGE  UNTIL  ZERO  HAS  BEEN  BRACKETED 

C 

C THIS  IS  AN  INTERNAL  SUBROUTINE  TO  EVALUATE  A  REDUCED  POLYNOMIAL 

C 
400   ICNT=N*2 
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P1=0.0 

401  Pl=Pl»Xl»A(ICNT-l) 
ICNT=ICNT-1 
IFUCNT-J)402.401,401 

402  60  TO  (320,330) . IP 


320 
321 


325 

326 

330 
331 

332 


360 


370 


380 
15 


41 


43 


1*1 


20 


30 


lF(P0«Pl)3?5.360t321 

X0  =  Xl 

X1=10.0«X1 

GO  TO  400 
■-ANSWER  LIES  BETWEEN  XO  AND  XI 

IP  =  2 

X2  =  X1 

Xl= (X0»X2) »0.5 

IF<ABS(X2-X0)-ABS(X1»TOLX) ) 360, 360 t 400 

IF(P0*PD331,360.332 

X2  =  X1 

GO  TO  326 

X0  =  X1 

GO  TO  326 
■-NOW  DIVIDE  OUT  THE  ZERO 

ROOTR( J-l ) =xl»DOM 

ROOTI (J-l ) =0.0 

I=N 

A(I)=A(I) *Xl*A(I*l) 

1  =  1-1 

!F(I-J)380.370,370 
J  =  JM 

CONTINUE 

DO  RO  J=J,N,2 

ITRY=0 

IF(J.LT.N)GO  TO 

P=A(N) 

0=A(N-1 ) 

GO  TO  50 

IF (NORD.GT.n)GO 


41 


TO 


43 

FOR 


STARTING  POINT 


TRY  PREVIOUS  ROOT 
P=PJ»R00TR (J) 

0=OJ«(ROOTR(  J)«»ROOTR(J>**OOTI  (J)  •ROOT  I  (J)  ) 
GO  TO  141 
P  =  0.0 
0=ITRY 

IF(ITRY,GT.l)GO  TO  47 
JTRY=ITRY*1 
DO  451  ITE»=1»ITMAX 
B0=0.0 
B1=0.0 
C0=0.0 
C1=0.0 
C2=0.0 

.-COMPUTE    COLUMN    OF    B'S    AND    C»S 
I=N*1 

OD=A(I)-P»B1-Q»BO 
B0  =  B1 
B1=DD 
1  =  1-1 

IF(I-J»2)40,40,30 
DD=B1-P*C2-Q«C1 
C0  =  C1 
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C1=C2 

C2=DD 

GO  TO  20 

C SPECIAL  FORM  FOR  LAST  COEFFICIENT 

40    C2=-P«C1-Q»C0 

0D=C1»C1-C?«C0 

DP=P0«C1-B1*C0 

DO=Bl«Cl-Pn«C2 

IF(A8S(00) ,LE.ABS(TOL?«DP).OR.ABS(00) .LE.ABS(T0L2«00) ) GO  TO  43 

OP=DP/OD 

DO=DQ/DD 

P=P«OP 

0=0*00 

ERR=DP»DP*OQ«DQ 

CRIT=P»P*0«0 

IF(CRIT.LE.T0L2«ERR)G0  TO  43 

CRIT=ERR/CRIT 

IFdTER.GT.1.  ANO.CRIT.LT. T0L31G0  TO  50 
C 00  AT  LEAST  ITMIN  ITERATIONS 

IF<ITEP-ITMIN><.5,<.5.44 
C IF  ERROR  IS  DECREASING.  KEEP  GOING 

44  IF  (CRIT.GE.CRITO)GO  TO  46 
IF«CRIT.LT.TOLT)GO  TO  SO 

45  ERRO=ERR 
CRITO=CRIT 

451   CONTINUE 
C DROP  OUT  WHEN  ERR  STARTS  TO  INCREASE  OR  ITER  IS  GREATER  THAN  ITMAX 

46  IF (CRITO-TOL)SO,50.43 

C IF  ACCURACY  IS  INSUFFICIENT.  GIVE  E3R0R  OUTPUT  AMD  P^OCE^D 

47  IR=3 
WRITE(5t5000)EDRO.ERR»CRlTO.CRIT 

SOOO  FORMAT (IP, 10G1?.4) 
SO    PP=-P/2.0 

DISC=PP«PP-0 

IF<DISC)52.5<..S4 
52    OISC=S0RT (-DISC) 

ROOTR ( J-l ) =PP«nOM 

ROOTR( J>  =Pd«oqm 

ROOTI (J-l ) =-DISC°DOM 

ROOTI (J) =DISC«DOM 

GO  TO  56 
54    DISC=SOPT (DISC) 

ROOTR (J-l )=(PP-DTSC)»DOM 

ROOTR (J) - (PP'DISO »00M 

ROOTI (J-l ) =0.0 

ROOTI (J)=0.0 

C SYNTHETIC  DIVISION  WITH  QUADRATIC  FACTORS 

56    I=N 

JL  =  J 

A(N)=A(N)-P 
60    1=1-1 

IF(I-J)80,70.70 
70    A(I)=A(I)-o»A(I»l)-Q<»A(T«2) 

GO  TO  60 
80    CONTINUE 

IF( JL-N)90. 100*100 
90    ROOTR(N)=-A(N) 

ROOTI (N) =0.0 
100   CONTINUE 
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r 


200  IER=IR 
RETURN 
END 
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APPENDIX  C  -  PROGRAM  USER'S  MANUAL 

This  appendix  will  discuss  the  calculations  performed  in 
the  various  subroutines,  the  output  format  options,  and  the 
input  data  deck. 

DESCRIPTION  OF  PROGRAM 

MAIN  -  The  Main  program  does  some  simple  calculations, 
contains  formats  for  I/O,  determines  from  the  I/O  options 
specified  what  calculations  are  to  be  performed,  and  calls 
the  appropriate  subroutines  to  perform  the  calculations. 

SUBROUTINE  UGLYDK  -  This  subroutine  takes  the  five 
inputed  offsets  and  uses  a  spline-cubic  curve  fitting  routine 
to  interpolate  ten  offsets  in  equal  increments  of  T'/IO  from 
the  base  line  to  T'.   The  subroutine  is  called  once  for  each 
set  of  inputed  offsets.   The  subroutine  requires  solution  to 
a  set  of  simultaneous  linear  equations  which  is  done  by  the 
IBM  SSP/SIMQ  subroutine  which  should  be  available  to  most 
users . 

SUBROUTINE  SPOO  -  This  subroutine  performs  the  calculations 
necessary  to  find  the  form  coefficient  k.   For  values  above 
k  =  0.19  equation  (6a)  is  used  to  determine  k.   This  subroutine 
performs  the  mapping  necessary  to  determine  the  regression 
analysis  form  coefficients.   This  subroutine  is  called  once 
for  each  set  of  inputed  offsets.   The  subroutine  requires  the 
roots  of  a  fifth  degree  polynomial  which  is  accomplished  in 
the  BAIRS  subroutine. 
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SUBROUTINE  SLAMPR  -  This  subroutine  takes  the  relative 
motions,  relative  velocities,  form  coefficients,  significant 
wave  height,  and  drafts  and  calculates  the  following: 

1.  number  of  slams  at  each  half-station, 

2.  longitudinal  extent  of  slam, 

3.  probability  of  slam  at  each  half-station, 

4.  most  probable  extreme  pressure  at  each  half -station, 

5.  extreme  pressure  for  design  at  each  half-station, 

6.  and  the  most  probable  and  extreme  velocities  of 
impact  at  each  half-station. 

This  subroutine  uses  the  various  equations  developed  in  Section 
IV  to  determine  the  above  information. 

SUBROUTINE  RELMV  -  This  subroutine  takes  the  non-dimen- 
sional relative  motions  and  relative  velocities  inputed  at 
station  1  and  4  from  the  tables  and  uses  linear  interpolation 
to  find  the  values  for  the  17  half-stations.   Seventeen  half- 
stations  were  considered  since  slam  damage  has  been  reported 
as  far  as  40%  aft  of  the  forward  perpendicular.   Linear  inter- 
polation is  valid  within  the  linear  theory  used  in  seakeeping 
programs . 

SUBROUTINE  FORCES  -  This  subroutine  performs  the  temporal 
integration  of  pressure  determined  by  SLAMPR  and  PRESS  sub- 
routines.  The  output  is  the  pressure  time  history  for  each 
half-station  in  the  slam  region. 
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SUBROUTINE  PRESS  -  This  subroutine  takes  the  pressure  at 
the  keel  from  SLAMPR  and  calculates  the  girth  distribution  of 
vertical  pressure  assuming  linear  variation  from  the  flat  to 
T'  and  resolves  the  vertical  component  at  each  elevation  using 
the  hull  offsets  to  calculate  the  angle  3  that  the  normal 
makes  with  the  vertical  (Figures  8a  &  8b) . 

SUBROUTINE  KEFFEC  -  This  subroutine  is  not  used  in  any 
of  the  calculations  but  determines  an  average  hull  form  coeff- 
icient.  This  value  is  a  measure  of  a  ship's  overall  suscep- 
tibility to  slam  forces.   For  two  similar  ships  the  one  with 
the  larger  average  k  will  have  larger  slam  forces  and  impulses. 
The  average  is  done  over  the  forward  eleven  half-stations 
since  slamming  rarely  occurs  aft  of  this  point. 

SUBROUTINE  APPROX  -  This  subroutine  uses  Simpson's 
integration  to  find  the  slam  impulse  by  integrating  the 
pressure  time  history  from  FORCES  at  each  half -station.   This 
subroutine  calls  the  IBM  SSP/QSF  subroutine  to  perform  the 
integrations.   Once  the  impulse  is  known,  the  amplitude  of 
an  equivalent  half-sine  pulse  with  equal  duration  can  be 
found.   This  sine  approximation  is  determined  to  input  into 
the  S.H.V.R.S.  dynamic  analysis  program. 

SUBROUTINE  MAXFOR  -  This  subroutine  fits  a  parabola 
through  the  three  highest  force  points  at  each  half-station  and 
the  first  derivative  is  set  equal  to  zero  to  find  the  maximum 
force  at  each  half-station. 
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SUBROUTINE  PLOT  -  This  routine  uses  the  M.I.T.  Joint 
Computer  Facility  PICTR  plotting  routine  to  plot  the  force 
time  history  for  each  half-station  on  one  graph.   Figures  9 
to  12  are  computer  plots  using  this  routine.   This  is  the  only 
subroutine  that  would  require  changing  if  running  the  program 
on  another  computer. 

OUTPUT  OPTIONS 

There  are  several  output  options  available  from  the 
program  depending  on  the  information  desired.   Options  are 
specified  in  the  data  deck  by  assigning  integer  values  to 
integer  variables.   The  following  options  are  available: 

1.  101  =  1   This  option  prints  information  from  the  SPOO 
subroutine  to  aid  in  determining  problems  in  this 
routine.   See  Figure  22  for  the  information  printed. 
101  =  0      Delete  the  above  output.   This  is  the 
normal  option. 

2.  103  =  1   This  option  allows  a  plot  of  the  pressure 
time  history.   See  Figures  13  to  16  for  typical  plots. 

103  =  0      Delete  the  plots. 

3.  104  =  1   This  outputs  a  table  of  inputed  offsets  and 
the  ten  calculated  offsets  for  each  half-station.   It 
is  desirable  to  specify  this  option  the  first  time  a 
set  of  lines  are  read  so  as  to  check  the  data.   See 
Figure  2  3  for  an  example  output. 

104  =  0      Delete  table  of  offsets. 
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4.  106  =  1   This  option  specifies  that  the  calculations 
and  output  values  will  use  most  probable  extreme 
pressure  values. 

106  =  0      This  option  specifies  that  the  program  uses 
the  extreme  pressure  for  design. 

5*  I07  =  !   The  force  time  history  is  outputed  in 
tabular  format.   See  Figure  24  for  an  example  of 
the  output  information. 

107  =  0   Delete  table  of  forces. 

6.  1010  =  0  Write  table  of  information  for  each  half- 
station  in  the  slam  region.   See  Figure  25  for  the 
information  available  in  this  option. 
1010  =  1  Write  a  single  table  of  information  for 
all  the  stations  in  the  slam  region.   Tables  7  to  10 
were  prepared  using  this  option. 
1010  =  2   Delete  both  of  the  above  tables. 

Any  combination  of  the  above  options  may  be  specified  to  get 
the  information  desired.   Since  the  main  cost  of  running  the 
program  is  the  output  cost,  it  was  felt  that  being  able  to 
eliminate  unnecessary  output  information  was  necessary. 
Several  input  options  are  available: 

1.  105  =  1   This  allows  the  inputing  of  all  the  relative 
motions  and  velocities  from  seakeeping  results  if 
available. 

105  =  0   Input  relative  motions  and  velocities  for 
stations  1  &  4  only.   (Non-dimensional  values  from 

Tables  [15]). 
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2.  Time  of  storm  may  be  specified  if  desired.   If 
TSTORM  is  specified,  it  will  be  used  in  the  calcul- 
ations.  If  set  equal  to  zero,  the  TSTORM  values  will 
be  obtained  from  Figure  4. 

3.  The  traveling  velocity  may  be  specified  if  desired. 
The  value  specified  should  be  for  a  520  foot  ship 
and  the  program  will  Froude  scale  this  figure.   If 
the  traveling  velocity  (VAFT)  is  set  equal  to  zero, 
the  traveling  velocity  will  be  calculated  by  the  pro- 
gram as  discussed  in  Section  IV. 

4 .  The  Froude  number  or  the  speed  in  knots  can  be 
specified.   Both  need  not  be  inputed. 

DATA  DECK 

In  this  section  the  input  format  for  the  data  deck  will 
be  discussed. 

CARD  SET  #1 

CARD  1  -  Format  (3  I  5)  -  101,  102,  104 

The  values  for  101  and  104  are  discussed  above. 

102  is  always  set  equal  to  1  and  indicates  that 

ships  lines  follow. 
CARD  2  -  Format  (20  A  4)  -  Title  and  information  on  ship. 

This  card  contains  any  information  desired  by  the 

user  to  designate  the  ship  being  run.   This  data  is 

printed  in  the  output  as  a  heading  only  and  is  not 

used  by  the  program. 
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CARD  3  -  Format  (5  F  10.5)  -  BB,  HH,  ALBP ,  ALWL,  STASP 
BB  =  Half  Beam 
HH  =  Design  Draft 

ALBP  =  Length  Between  Perpendiculars 
ALWL  =  Length  of  Load  Water line 
STASP  =  Station  Spacing  (not  half-stations) 

CARD  SET  #2 

This  card  set  contains  the  ships  lines  necessary  for 
slam  calculations. 

CARD  1  -  Format  (5  F  10.5) 

This  card  contains  the  waterlines  at  which  the 
offsets  were  taken.   The  first  value  must  be 
the  0.0  foot  waterline. 

CARD  2  -  Format  (5  F  10.5) 

The  offsets  for  the  five  waterlines  in  the  same 
order  as  CARD  1  for  the  first  half-station. 
CARD  2  is  repeated  for  as  many  half-stations  as 
desired  up  to  a  total  of  17  half-stations. 

CARD  3  -  BLANK  CARD 

This  card  indicates  that  the  last  set  of  offsets 
has  been  inputed. 

CARD  SET  #3 

This  set  contains  only  one  card.   If  no  plots  are  desired 
(103  =  0) ,  this  card  is  left  blank. 
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Format  (I  5,  4  F  10.5)  -  IPLT,  XSCL 

IPLT  =  +1  -  This  specifies  that  the  plot  routine 

scales  the  plots.   (AUTOSCALE) 

IPLT  =  -2  -  This  specifies  to  the  plot  routine  that 

the  values  in  XSCL  be  used  to  scale  the  plots. 

XSCL(l)  =  0.0  =  Time  slam  starts 

XSCL (2)  =  End  of  time  axis 

XSCL (3)  =0.0 

XSCL (4)  =  Maximum  force  on  vertical  axis. 
If  IPLT  =  +1,  XSCL  need  not  be  specified.   This  same 
data  will  be  used  to  scale  all  plots  until  a  new  set  of  lines 
are  read. 

CARD  SET  #4 

Card  set  #4  contains  the  data  necessary  for  a  particular 
set  of  operating  conditions.   This  set  is  repeated  for  as  many 
operating  conditions  as  desired. 

CARD  1  -  Format  (6  I  5)  -  109,  103,  105,  106,  107,  1010 

109  =  1   always  =  1  -  indicates  a  new  set  of  operating 
conditions  follow.   The  other  integer  variables 
were  discussed  above. 
CARD  2  -  Format  (20  A  4) 

This  card  contains  information  which  is  printed  on 
the  output  to  distinguish  this  operating  condition 
from  the  others. 
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CARD  3  -  Format  (8  F  10.5)  -  FROUDE,  SPEED,  SEA,  HP,  HA, 

CB,  VAFT,  VTHR 

FROUDE  =  Froude  number 

SPEED  =  Ship  speed  in  knots 

SEA  =  Significant  wave  height  (feet) 

HF  =  Draft  Forward  (feet) 

HA  =  Draft  Aft  (feet) 

CB  =  Block  coefficient  of  underwater  portion  of 
hull  for  this  draft. 

VAFT  =  Traveling  velocity  for  520  foot  ship  (ft/sec) 

VTHR  =  Threshold  velocity  for  520  foot  ship  (ft/sec) 
CARD  4  -  (F  10.5)  -  TSTORM 

TSTORM  =  Duration  of  storm  in  hours 
CARD  5  -  IF  105  =  0,  Format  (4  F  10.5)-  RM(2),  RV(2), 

RM ( 8 ) ,  RV ( 8 ) 

RM(2)  =  relative  motion  from  tables  at  station  1 

RV(2)  =  relative  velocity  from  tables  at  station  1 

RM(8)  =  relative  motion  from  tables  at  station  4 

RV(8)  =  relative  velocity  from  tables  at  station  4 

IF  105  =  1   Format  (8  F  10.5)  for  each  card  -  use 
as  many  cards  as  necessary  to  specify  the  non-dimen- 
sional relative  motions  and  velocities  for  as  many 
half-stations  as  offsets  read  in  Card  Set  #1.   The 
values  are  specified  in  pairs  beginning  at  half- 
station  1  (i.e.,  RM(1),  RV(1),  RM(2),  RV(2),  ) 
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CARD  6  -  IF  103  =0—  Delete  Card 

IF  103  =1  —  FIRST  40  spaces  contain  X-axis 
label 

NEXT  40  spaces  contain  Y-axis 
label 
This  card  is  used  by  plotting  subroutine. 

At  the  end  of  the  operating  conditions  for  this  ship  a 
blank  card  is  inserted  to  indicate  the  end  of  the  calculations 
for  this  set  of  lines.   The  above  card  sets  1  to  4  can  be 
repeated  for  as  many  ships  as  desired.   After  the  last  ship 
has  been  read,  another  blank  card  is  inserted  to  indicate  the 
end  of  the  data  deck.   Therefore  the  last  two  cards  in  the  deck 
will  be  blank  cards  indicating  the  end  of  the  operating  condi- 
tions and  that  the  last  ship  has  been  read. 
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FIGURE    23   -    LIST  OF    INPUTED    AND    CALCULATED   OFFSETS    (IQ4=1) 
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* 


* 

*•&*$  *  156-7  » 

J».070  *  43  3.  P  * 


■7*^?  *  f36«5 


*  <"•»*/  *  ^-36.5  * 


* 


7!  •  7  <}  7  «  ^p,q.q 


* 


*  ^  •  1  *  ?  *  1  7  fc  3  •  3  « 

*  Pi  •  1 1  /  *  ^  1 3  •  A  « 
£•127  *  7P4.1  * 

*  71  •  1  3  0  *  A  8  7  •  6  * 
»                    Pi  •  1  **  ?  w  2  P  7  •  2  * 

*  7.  •  1 5  7  *  3  3  •  8  * 
v  '"•167  *  7  •  I?  * 
»                       7  •  1  7  \  m  0  •  C  * 

*  P>  •  1  8  1  *  7i  .  7  « 
?-  •  *  9 1  *  (? .  7  « 
^•2P)1  *  I?. .7  « 

*  3 • 2 1 3  *  7.7  * 
»                        »  •  2  ?  1  *  p. .  7  * 

*  7-  •  2  3 1  *  7 .  7  „ 
0*2^1  *  0.?  * 


FIGURE    24   -   FORCE    TIME    HISTORY      (  IQ7  z  1  ) 
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*  *****************************  ********** 

*  * 

*  HAl  F     STATION!  h  * 

*  * 
**************************************** 
rXTPEMF     PRESSURE     =                         ?4/».     PSI 

«OST     Pn'QRARLF     oPrpsUPF     =  iff.     PSI 

rOPK<     C^FF.      (K)     =  C*.^7434 

VOM-DIVEMSIONAL     PFLATIVE     MOTTON     =        Pl.jflglR 
vCNJ-OIwEMSIONAL    PFLATIVF    VELOCITY    =       H.tf8K0 
vnv»pp    nr     si.ANS    DURING    EXTPE^P     STORM    =  k7XH. 

nPO^AEJILlTY     OF     SLAM     IMPACT     =       Pl.?7?<»? 
rXT'-Tir     RELATIVE     VELOCITY    =  57.3    F7/SFC 

V.OCT     PPORABl„F     RELATIVE     VELECTTY     =  47.?    FT/SEC 

MAXIMUM    VALUE     Or    SLAM    FORCE    *»  3543#0    MPS 

TMPH|_cr     VALUE     =  183.     K  T  P  -  S  E  t 

AMPLITLIDF    Oc     SINF     APPROXIMATION     =  P37P.7    KIPS 

nURATI^N    Or     Si  A^     IMPACT    FO^CF     =        C-.t?'/9?    SEC 
TJVE     Si   AM     STARTS     =        w,t^R86     Sf C 


FIGURE    25 -INFORMATION    PRINTED    AT   EACH    HALF-STATION 
(IO10  z  <t» 
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APPENDIX  D  -  EXAMPLE  PROGRAM 

An  example  program  data  deck  and  sample  output  are 
included  to  aid  in  using  this  program.   MARINER  case  2  is 
the  example  run  which  is  the  fully  loaded  case  at  15.5  knots 
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//  XEO  SLAMCW 


0 

1 

0 

M  AD  T  kJr~  ri 

3R.0 

™  AK  i  rsj 

29.75 

r  « 

528.0 

537.13 

26. 

0.00 

1.00 

2.00 

3.00 

4.00 

0.5 

1.5 

1.75 

2.05 

2.3 

1.00 

3.00 

3.5 

4.1 

4.6 

1.0 

3.45 

4.3 

5.15 

5.6 

1.0 

3.9 

5.1 

6.2 

6.6 

1.5 

4.95 

6.T 

7.5 

8.05 

?.6 

6.0 

7.S 

8.8 

9.5 

4.0 

7.5 

9.2 

10.6 

11.6 

5.4 

9.0 

10.9 

12.4 

13.5 

6.95 

11.0 

13.1 

14.8 

1A.0 

8.00 

12.00 

15.00 

17. 

18.2 

8.58 

15. 

17.6 

19.7 

21.0 

9.17 

18.00 

20.2 

22.3 

24.08 

13.6 

21. 

23.10 

24.9 

26.5 

-2 

0, 

.0 

0.5 

0. 

,0 

5000 

• 

1 

1 

0 

0 

0 

1 

2 

MARINER    -    F 

=  0. 

20       S=0, 

,047       C3=0, 

,613 

15.5 

0.047 

27.0 

27.0                  0.613 

0.0 

0.0280 

0.095 

0.0191 

0.069 

2 

MARINER 

TIME 

(SFO 

FO 

BLANK 


fully  loaofd 

0.0  12.0 


(KIPS) 

BLANK 
BLANK 


//    END 


FIGURE    26a  -    SAMPLE     PROGRAM   -   DATA     DECK 
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THIS  PROGRAM  CALCULATE*  THE  SLAM  FORCE  DISTRIBUTION  VS  TIME 
ON  THI-  LOWEK  1/10  OF  THE  SHIP  FROM  STATION  8  FORWARD. 


MAkTN<"R AS  BUILT 

uALp  BfAM  .    ^8.tlPirfiPl  DFSIGN  DRAFT  -    20.75*00 

i'BP  ■      S?8.*0  LENGTH  LOAD  WATER  LINE  -      537.13 

Distance  betwefn  stations  -  pa.ap 

(ALL  DISTANCES  IN  FEETi 


AVtRAGE  HULL  SHAFE  FACTOR  -   0.11909 


FIGURE    26b   -    SAMPLE     PROGRAM     -   OUTPUT 
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